Preprint typeset in .THEP style - HYPER VERSION 

CCTP-2010-6 
ITP-UU-10/19 
SPIN-10/17 

Langevin difFusion of heavy quarks in 
non-conformal holographic backgrounds 



Umut Gursoy\ Elias Kiritsis^"^, Liuba Mazzanti^ and Francesco Nitti^ 

^Institute for Theoretical Physics, Utrecht University; Leuvenlaan 4 , 3584 CE 
Utrecht, The Netherlands. 

^ Crete Center for Theoretical Physics, 

Department of Physics, University of Crete 71003 Heraklion, Creece 

^ Departamento de Ftsica de Particulas, Universidade de Santiago de Compostela 
and Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E- 157 82, Santiago de Compostela, Spain 

^APC, Universite Paris 7, Bdtiment Condorcet, F-75205, Paris Cedex 13, 
France, (UMR du CNRS 7164). 

Abstract: The Langevin diffusion process of a relativistic heavy quark in a non- 
conformal holographic setup is analyzed. The bulk geometry is a general, five- 
dimensional asymptotically AdS black hole. The heavy quark is described by a 
trailing string attached to a flavor brane, moving at constant velocity. From the 
equations describing linearized fluctuations of the string world-sheet, the correlation 
functions defining a generalized Langevin process are constructed via the AdS /CYT 
prescription. In the local limit, analytic expressions for the Langevin diffusion and 
friction coefficients are obtained in terms of the bulk string metric. Modified Einstein 
relations between these quantities are also derived. The spectral densities associated 
to the Langevin correlators are analyzed, and simple analytic expressions are ob- 
tained in the small and large frequency limits. Finally, a numerical analysis of the 
jet-quenching parameter, and a comparison to RHIC phenomenology are performed 
in the case of Improved Holographic QCD. It is shown that the jet-quenching pa- 
rameter is not enough to describe energy loss of very energetic charm quarks and the 
full Langevin correlators are needed. 
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1. Introduction and results 

RHIC experiments of heavy-ion collisions and related data on the deconfined phase 
of QCD, [1] have provided a window for string theory techniques to meet the real 
world. The context is strong coupling dynamics near and above the deconfining 
transition in QCD. String theory via the AdS/CFT correspondence has provided a 
framework in order to understand strong coupling dynamics in the deconfined phase 
including the calculation of transport coefficients. Recent reviews on the progress in 
this direction are [2, 3, 4, 5, 6]. 

Observables of particular importance are associated to heavy quarks. Heavy 
quarks may be produced in the Quark-Gluon Plasma (QGP) of the RHIC fireball 
and are then travelling to the detectors while moving through the dense QGP. They 
can be tagged reasonably well and are therefore valuable probes of the dynamics in 
the plasma and in particular for the mechanism of energy loss. 

A single heavy quark can be modeled in string theory by an open string. Its 
end-point is representing the heavy quark while the string is trailing behind as the 
quark moves. The large mass limit is important in order to neglect the non-trivial 
flavor dynamics associated with light quarks (although with improved techniques the 
light quarks may also eventually be addressed reliably in the holographic context). 
As quarks are associated with strings ending on flavor branes, a heavy quark ends on 
a brane that is stretching in the UV part of the bulk geometry. The motion of such a 
string, and the associated force acting on the quark from the thermal medium, have 
been studied in detail with several complementary methods, [7, 8, 9]. In the simplest 
setup, the UV endpoint of a fundamental string is forced to move with constant 
velocity v along a spacial direction. The equations of motion for the full string are 
solved and the radial proflle of the trailing string is found as it moves in a bulk black- 
hole background representing the deconflned heat bath. The energy absorbed by the 
string is calculated and the drag force of the string is obtained. The picture remains 
roughly valid, while details change when conformal invariance is broken, [10, 11]. 

An important improvement in this picture consists of the study of the stochastic 
nature of this system in analogy with the dynamics of heavy particles in a heat 
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bath giving rise to Brownian motion. This involves a diffusive process, that was first 
considered in a holographic setting in [12], by using the Schwinger-Keldysh formalism 
adapted to AdS/CFT in [13]. 

Subsequently, a study of the (quantum) fluctuations of the trailing string, [14, 15] 
provided the information on the momentum broadening of a heavy quark as it moves 
in the plasma. The stochastic motion was formulated as a Langevin process, [16, 17] 
associated with the correlators of the fluctuations of the string. 

Many heavy quarks in experiments are relativistic. Therefore it is necessary 
to study the associated relativistic Langevin evolution of the trailing string, a feat 
accomplished in the A/" = 4 case in [18]. The same type of Langevin process was 
studied in [19] for the case of an accelerating quark in the vacuum (rather than in a 
deconfined plasma), by analyzing the fluctuations of a trailing string in AdS with a 
non- uniformly moving endpoint. 

On the experimental front, there have been several results from the RHIC exper- 
iments, [21]-[25]. The experimental signatures are currently summarized by the 
spectra that originate in the semileptonic decays of charmed and bottom hadrons. 
From these spectra a modification factor i?^^ and an elliptic flow coefficient ^i^s 
extracted. They capture the effects of the medium to the propagation of the heavy 
quarks. The data exhibit a substantial elliptic flow, up to ~ 10%, and a high-pr 
suppression down to i?^^ ~ 0.25. These values are comparable to light hadrons. 
Radiative energy loss models based on pQCD, [26] do not seem to explain well the 
experimental data, [27]. Elastic scattering energy loss plus non-perturbative interac- 
tions can on the other hand accommodate the data, [27]. 

In particular the Langevin approach has been applied to the study of the heavy 
quark energy loss by several groups, and the related physics is summarized in the re- 
view [27]. The Langevin evolution used was relativistic and with symmetric diffusion 
coefficients. As there was no microscopic model to provide the proper fluctuation- 
dissipation relation, the Einstein equations used vary, and in all examples it was 
assumed that the equilibrium distribution is the Jiittner-Boltzmann distribution. 
Moreover various combinations of friction forces were used, resonance models, pQCD, 
N=4 AdS/CFT and combinations. A further recent analysis was performed in [28] 
with similar conclusions. The associated relativistic and isotropic Langevin systems 
used have been introduced in the mathematical physics literature rather recently, 
[29] (see [30] for a review). 

The purpose of the present paper is to study further the relativistic Langevin 
evolution of a heavy quark using holographic techniques in a general context, and 
going beyond conformally invariant backgrounds characterized solely by AdS geome- 
tries^ . In this work we will study a large class of non-conformal backgrounds captured 
by Einstein-dilaton gravity with a dilaton potential in 5 dimensions. In a series of 

"'^A particular example in this class was recently considered in [20], which studied the Langevin 
process in the non-conformal Af = 2* background. 
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recent works, such backgrounds were analyzed both quahtatively and quantitatively 
and have provided a rich variety of holographic bulk dynamics. In particular, for a 
selected class of scalar potentials, they mimic the behavior of large- Yang Mills, 
[32]-[38]. This match can be quantitative, [38], agreeing very well both at zero and 
finite temperature with recent high-precision lattice data, [39]. On the other hand, 
the analysis we carry out in the present paper is general, as it applies to any asymp- 
totically AdS background. 

We therefore consider a fundamental string whose end-point lies in the UV region 
of a bulk black-hole background of a non-conformal holographic model. The string 
end-point is forced to move with velocity v. Solving the Nambu-Goto equations of 
motion, the classical profile of the trailing string can be found. The string stretches 
inside the bulk until it becomes completely horizontal at some value of radial coordi- 
nate Ts, given by /(r^) = v"^ where /(r) is the blackness function of the background. 
When the quark is moving slowly, as w — )■ 0, the point approaches the bulk black 
hole horizon. 

The induced metric on the string world-sheet has the form of a two-dimensional 
black-hole metric with a horizon at r = as first observed in [15].^ This black- 
hole is an important ingredient of the dynamics of the system. In particular it is 
crucial in the calculation of the thermal correlators using the Schwinger-Keldysh 
formalism, as well as for the fluctuation-dissipation relation. The world-sheet black 
hole has an associated Hawking temperature Tg that depends on several parameters: 
the background temperature T, the zero-temperature bulk scale and the quark 
velocity v. It coincides with the temperature T of the heat bath only in the non- 
relativistic limit. In the conformal case, one has Tg = Tsconf = ^(1 ~ v"^)^ ^ T. 
The numerical analysis performed in Section 6 shows that, in our non-conformal 
holographic model < Tg^conf < T. The equality Ts = Ts^conf, in the first relation is 
attained, for arbitrary v, in the high T limit, as shown in figure 1 in section 6 and 
also in the ultra- relativistic limit, w — )■ 1. 

We next consider small fluctuations around the classical string profile. They sat- 
isfy second-order radial equations that are related to the associated thermal correla- 
tors by the holographic prescription. It should be emphasized that such correlators 
are thermal with temperature and not the temperature T of the heat bath^. More- 
over, they satisfy the fluctuation-dissipation relation associated with the emergent 

^This is a generic effect on strings and D-branes embedded in black-hole/black-brane back- 
grounds. It was first observed in [40] where it was used to propose that a different speed of hght 
is relevant for such branes. It is impUcit or expUcit in many holographic computations using probe 
flavor branes, [41, 42] and strings [15]. 

^This appears as an integration constant in the background geometry and corresponds to the 
dynamically generated energy scale Kqcd in the dual field theory. 

"^More precisely put, the Langevin correlators that are obtained from the string fluctuations by 
the holographic prescription obey a modified Einstein relation with temperature rather than T. 
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temperature Tg. The fact that the string fluctuations see a modified temperature 
crucially affects the Einstein relation between the diffusion constants. 

At the quadratic level of fluctuations, a relativistic Langevin diffusion equation 
is obtained using the AdS /CFT prescription. The diffusion constants and friction 
coefficients are calculated analytically in terms of the bulk metric, for general non- 
conformal backgrounds. The Einstein relation is now modified, due to the fact that 
the temperature is modified. Another important property is that the diffusion con- 
stants perpendicular and longitudinal to the motion (denoted by k± and are 
different, a fact that was already observed in the conformal relativistic case, [f8]. 
This is persisting here, and we are able to show, for general non-conformal back- 
grounds, that K\\ > K±, namely the longitudinal diffusion constant always dominates 
the transverse one. Furthermore, both the diffusion constants and the friction term 
are momentum dependent, as expected. This is in contrast to the conformal case. 

The properties of this relativistic Langevin evolution differ substantially from 
rotationally invariant equations that have been introduced recently in mathematical 
physics [29, 30, 28]. In particular, here the evolution is not rotationally symmetric, 
and the Einstein relation is different, because the fluctuation-dissipation relation 
is different. This implies that the equilibrium configuration is not the standard 
rotationally-invariant Jiittner-Boltzmann distribution. 

The processes we discuss here are connected to a general property of a class 
of statistical systems, and provide a concrete and rather general solvable example 
thereof. Such systems, when in contact with a heat bath of temperature T, if stirred 
gently and continuously, end up in a stationary state that is thermal but with a 
temperature Tg, different from that of the heat bath. They satisfy a fluctuation- 
dissipation relation involving the new temperature, [43] . This phenomenon has been 
expected to occur in general, in situations with slow dynamics. These include in 
particular glassy systems, as well as systems that are very gently stirred by external 
agents and reach stationarity. Here, we have a system that is stationary but strongly 
driven by the external source (an electric field on the flavor brane, that keeps the 
velocity large and constant as the quark is moving through the plasma). Typically, 
Tg > T, [43] but in our case things are different. In the conformal case Tg = Tg^conf = 
T{l-v^)i < T and this seems to persist in several other cases, as already mentioned 
before. 

The problem we are solving can be also cast in a different light. The follow- 
ing question has been asked since the beginning of the 20th century: what are the 
Lorentz-transformation properties of temperature and the associated stationary dis- 
tribution? This question is still considered open, [44] with several conflicting answers. 
Our setup can be reconsidered as follows: the heavy quark is a "thermometer" mov- 
ing inside a heat bath of temperature T. The way it measures temperature is via the 
fluctuation-dissipation theorem as argued in generality in [43]. Therefore the tem- 
perature that it measures as it moves in the thermal medium is Tg, which is velocity- 
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dependent. This velocity dependence is simple in the conformal case, Tg = T{l — v'^)^ 
but is rather complicated in the non-conformal case and is therefore system- (and 
possibly thermometer-) dependent. 

The Langevin correlators must be renormalized as they are divergent near the 
AdS boundary of the string world-sheet. We show that the only UV divergence 
is subtracted by a counterterm that renormalizes the (heavy) quark mass. The 
associated scheme dependence affects the real part of the correlators only. 

The local Langevin equation arises when looking at the large-time limit of the 
fluctuations of the heat bath, i.e. at the small frequency modes. On the other 
hand, the holographic computation gives access to the full frequency spectrum of 
the correlation functions driving the generalized Langevin dynamics. In this work 
we compute holographically the full Langevin correlators and the associated spectral 
densities. In particular, we obtain analytic expressions (in terms of the bulk metric 
and dilaton profiles) in the two opposite regimes of small and large frequencies uj 
(compared to an appropriate temperature scale). 

For small frequencies, an analytic expression for the spectral densities is obtained 
using the membrane paradigm [45], which allows to relate these quantities to the near- 
horizon values of the background functions. In the large-frequency regime, on the 
other hand, the spectral densities are obtained via a modified WKB method, similar 
to the one followed in [46] for bulk fluctuations in an AdS'-Schwarzschild background. 
The high-frequency behavior is different, depending on the mass of the probe quark. 
For finite mass, and for large u, the spectral densities grows linearly with uj, whereas 
in the limit when the quark mass becomes infinite this behavior changes to a cubic 
power-law. 

Going beyond the zero-frequency limit is necessary when the diffusion process 
happens on time scales comparable to, or smaller than the auto-correlation time of 
the fluctuation propagators. More specifically, since these are thermal correlators 
at the temperature Tg, the large-time approximation breaks down over time-scales 
shorter than T~^. This condition puts a temperature-dependent upper bound on the 
momenta of the heavy quark, above which the diffusion process cannot be described 
by a simple local Langevin equation with white noise. In this context, it is useful to 
have approximate expressions for valid for large frequencies (for examples, those we 
obtain with the WKB method) to model the behavior of the system in the regime 
where the local Langevin approximation breaks down and the dynamics becomes 
non-markovian (due to a non-trivial memory kernel) 

The results described above apply to any five- dimensional holographic model 
which admits asymptotically AdS black-hole solutions. On the other hand, it is 
interesting to perform a quantitative comparison between the diffusion constants cal- 
culated in concrete models, and characteristic observables in heavy ion experiments 
that can simply be connected to Langevin processes. In the context of heavy-ion 
physics, the transverse diffusion constant is directly related to the jet-quenching pa- 
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rameter q± = 2k±/v. The latter is a convenient quantity to describe the observed 
phenomenon of transverse momentum broadening of a heavy quark: this is the pro- 
cess by which the transverse momentum of the heavy quark probe^, initially equal 
to zero, undergoes a stochastic diffusion process such that after a time t it acquires 
a dispersion Ap^ = q±vt. 

In the final part of this work we perform a quantitative analysis of both the full 
Langevin correlators, and of the jet-quenching parameter, in a particular Einstein- 
dilaton model, namely Improved Holographic QCD [32, 33], which agree quite well 
both qualitatively and quantitatively with the zero- and finite-temperature Yang- 
Mills theory. In particular, we focus on the specific model which was put forward in 
[38], and displays a good quantitative match with the spectral and thermodynamic 
properties of lattice Yang Mills theory. 

The analysis is performed numerically, both with respect to the background 
metric and to the solution of the fluctuation equations. By a shooting technique, we 
determine the wave-functions describing the world-sheet fluctuations, and obeying 
the appropriate retarded boundary conditions. From the wave- functions, the holo- 
graphic prescription allows to determine the full Langevin retarded correlator, whose 
imaginary part gives the associated spectral density. 

Using the exact numerical evaluation we are able to test the different analytic 
results discussed above. In particular, we test the validity of the WKB result for 
large frequency, and in various regimes of quark mass and velocity. Unexpectedly, 
we find that the analytic WKB formulae not only capture the large frequency regime, 
but are a very good approximation to the correlators at almost all frequencies. 

The numerical evaluation of the diffusion constants may lead directly to a com- 
parison of the jet-quenching parameters between the holographic QCD model and 
data. We find that q± displays a mild momentum dependence for large quark mo- 
menta, which however differs from the one obtained holographically in the conformal 
case. As the temperature rises, q± increases significantly, approximately as ~ 
Interestingly, it is found that for temperatures above ~ 400 MeV, the local de- 
scription of the diffusive process breaks down for charm quarks with momenta above 
~ 5 — 10 GeV. This is because the process occurs on time scales shorter than l/Tg. 
This would imply that in order to describe heavy charm quark diffusion in the AL- 
ICE experiment, one would need the full generalized non-local Langevin equation, 
and the full frequency-dependent correlator, rather than just its low-frequency limit 
captured by qj_. This would constitute an interesting testing ground for holographic 
models, where the full correlators can be easily computed. We also estimate the 
energies at which the energy loss mechanism described here, is not any more the 
dominant one and radiation becomes the dominant mechanism. This is estimated by 

^Here, "transverse" refers to the initial quark trajectory, not to the direction of the coUiding 
beams 
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requiring that remains below the would-be position of the flavor brane, [15, 53]. 
We show that these limit do not substantially constrain this framework. 

A direct quantitative comparison of the results of this paper with data is ham- 
pered by the fact that quark degrees of freedom in the plasma are not included in our 
analysis. We are compensating (partly at least) for this using the "energy scheme" 
for comparison, however the recent results found in [54] (and reviewed in [55]) suggest 
that even in that case we may be underestimating the result. 

In summary the AdS/CFT calculation of the Langevin diffusion of heavy quarks 
has the following characteristics 

• The diffusion coefficients are asymmetric. The longitudinal diffusion coeffi- 
cients is always larger than the transverse one. They both become large with 
increasing 7. 

• The Langevin correlators satisfy a thermal fluctuation-dissipation relation with 
a temperature Tg that is typically smaller than then heat bath temperature T. 
In the conformal case Tg = [18]. The associated Einstein relations are 
non-standard, especially the longitudinal one. 

• The local (Markovian) Langevin diffusion breaks down at some energy scale. 
Beyond this scale, the full force correlators are needed. This breakdown is 
expected to be relevant at LHC energies for the charm. 

This paper is organized as follows. Section 2 describes in detail the Langevin 
equation for a relativistic heavy quark travelling in the quark-gluon plasma. In this 
section we also review how to describe the Langevin dynamics in the holographically 
dual geometry, in terms of fluctuations of trailing strings. Section 3 presents the 
necessary background for the holographic computation. In particular we present the 
holographic dual geometry of our non-conformal model, the relevant classical trailing 
string solution, and the corresponding linear fluctuations. It is in this section that 
we obtain the fluctuation equations in general non-conformal black hole space-times, 
whose solutions enter the construction of the Langevin propagators. 

Sections 4 and 5 contain our main results. In Section 4 we discuss the Langevin 
correlators and the associated spectral densities, first in full generality, then in the 
various limits of low- and high- frequency. In section 5 we specialize to the low- 
frequency modes, which compute the long-time behavior of the diffusion and friction 
coefficients of the local Langevin equation. We provide exact analytic expressions 
for these quantities, in terms of the background metric functions. We also discuss 
the non-relativistic and ultra-relativistic limits, and derive the modified Einstein re- 
lations. While the previous sections deal with a completely general holographic dual, 
in Section 6 we provide a numerical study of these results in a specific model, namely 
Improved Holographic QCD that was shown in [38] to provide a good quantitative 
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description of the static properties of pure Yang-Mills at zero and finite temperature. 
In particular, in this section we compute the jet-quenching parameter arising from 
this model, and discuss the results in light of RHIC data. 

Several technical details are left to the Appendices. In Appendix A we dis- 
cuss some subtleties in the definition of the propagator, related to boundary terms. 
Appendix B provides the details of the calculation of the diffusion constants; In Ap- 
pendix C we give a detailed discussion of the WKB method that we use to obtain the 
large-frequency limit of the spectral densities, as it is more involved than the con- 
formal case. In Appendix D we discuss the Langevin correlators in the conformal, 
A/" = 4 case. 

2. Langevin Equation for a relativist ic heavy quark 

In this section we review how the diffusion of a relativistic heavy quark through the 
plasma is described by a generalized Langevin equation. First we give the purely 
4D picture. Then, in subsection 4.2, we review the holographic description of the 
Langevin process that appeared in the previous literature, for the case of AdS- 
Schwarzschild black- holes. This will be extended to general asymptotically AdS 
geometries in Section 3. 

2.1 The Langevin equation in the boundary theory 

Consider a quark which, in a first approximation, experiences a uniform motion across 
the plasma, with constant velocity v. Due to the interactions with the strongly- 
coupled plasma, the actual trajectory of the quark is expected to resemble Brownian 
motion. To lowest order, the action for the external quark coupled to the plasma 
can be assumed, classically, to be of the form: 



where 5*0 is the free quark action, and J^{t) depends only on the plasma degrees of 
freedom, and plays the role of a driving force (the "drag" force). 

To obtain an equation for the quark trajectory one needs to trace over the plasma 
degrees of freedom. If the interaction energies are small compared with the quark 
kinetic energy (therefore for a very heavy quark, and/or for ultra- relativistic prop- 
agation speeds), tracing over the microscopic degrees of freedom of the plasma can 
be performed in the semiclassical approximation, and the quark motion can be de- 
scribed by a classical generalized Langevin equation for the position X*(t), of the 
form: 






% = 



1,2,3 




9 



Here, C^^{t) is a memory kernel, 6{t) is the Heaviside function and ^(t) is a Gaussian 
random variable with time-correlation: 



{e{tW{t')) = A^'it-t') (2.3) 

The functions A'^^{t) and C^^{t) are determined by the symmetrized and anti- 
symmetrized real-time correlation functions of the forces J^{t) over the statistical 
ensemble: 

C'Kt) = Gly^mit) ^ -^{[^{t),FK^), A^^t) = G%^{t) ^ -l({^^(t),^^(0)}). 

(2.4) 

The results (2.2) and (2.4) are very general, and do not require any particular 
assumption about the statistical ensemble that describes the medium (in particular, 
they do not require thermal equilibrium). One way to arrive at equation (2.2) is 
using the double time formalism and the Feynman- Vernon influence functional [48] . 
A clear and detailed presentation can be found in [49], chapter 18. 

The retarded and advanced Green's function are defined by: 

Gl{t) = e{t)CHt), Gtit) = -e{-t)CHt), (2.5) 
which lead to the relation 

CHt) = G%{t) - Gtit) (2.6) 

Notice that the kernel entering the first term on the right in equation (2.2) is the 
retarded Green's function, G%{t) = e{t)C'^{t). 

It is customary to introduce a spectral density p''-'{uj) as the Fourier transform of 
the anti-symmetrized (retarded) correlator, 

/ + 00 p + oo ij / l\ 

dup'^{uj)e-"'\ G%iu)= du' ^ \\ . (2.7) 

From equation (2.6) and the reality condition GA{t) = GR{—t), or in Fourier 
space, Ga{(^) = G*i^{uj), we can relate the spectral density to the imaginary part of 
the retarded correlator: 

p'\uj) = --lmG%{u) (2.8) 

Local limit. Suppose the time-correlation functions vanish for sufficiently large 
separation, i.e. for times much larger than a certain correlation time Tc- Then, in 
the limit t ^ Tc, equation (2.2) becomes a conventional local Langevin equation, 
with local friction and white noise stochastic term. Indeed, in this regime the noise 
correlator can be approximated by 

A'^{t-t')^ K,'^S{t-t'), t-t'->Tc. (2.9) 
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This equation defines the Langevin diffusion constants Similarly, for the friction 
term, we define the function Y^{t) by the relation: 

cHt) = j^iHt) (2.10) 

so that the friction term can be approximated, for large times, as: 

/OO / POO \ 

dTC'^iT)Xjit -r)-U dr^'^ir)] Xjit), t » r,. (2.11) 

In this regime, equation (2.2) becomes the local Langevin equation with white 
noise, 

^ + v''x,{t) = at), {cm'it')) = K^'s{t - 1% (2.12) 

with the self-diffusion and friction coefficients given by: 

k'' = lim G%^{u)- 7]'' = / drY^ir) = - lim (2.13) 

In the case of a system at equilibrium with a canonical ensemble at temperature 
T, one has the following relation between the Green's functions: 

Gsymi^) = -coth^lmGRiu), (2.14) 

which using equation (2.13) leads to the Einstein relation k,^^ = 2Tri^^ . For such a 
thermal ensemble, the real-time correlators decay exponentially with a scale set by 
the inverse temperature, therefore the typical correlation time is Tc ~ 1/T. 

Determining and studying the Langevin correlators (2.4) , and the diffusion 
constants (2.12) will be the main purpose of the rest of this paper. 

Next, we write down explicitly the classical part, 6SQ/6X{t) of the Langevin 
equation, in order to arrive at an equation describing momentum diffusion. We start 
with the kinetic action for a free relativistic quark. 



Sm^-M^l^r^'-^f (2,15) 

We choose the gauge r = X°, and obtain 

6So/6X\t) =dpi/dt , with , pi = MgXi{l-XiX'y^/'^. (2.16) 
Equation (2.12) becomes the Langevin equation for momentum diffusion: 

% = -V%{p>,+e{t), (2.17) 



where: 



vliP ) = IT^fW' ^(p ) = \/1+pVMI (2.18) 
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Linearized Langevin equations. For a generic quark trajectory, the Langevin 
equation (2.17) is non-linear, due to tlie p-dependence implicit in 77^. To put it in 
a form which allows for the holographic treatment in terms of the trailing string 
fluctuations, it is convenient to derive from equation (2.17) a linearized Langevin 
equation for the fluctuations in the position around a trajectory with uniform ve- 
locity, X(t) = vt + 5X. To this end, we separate the longitudinal and transverse 
components of the velocity fluctuations: 

X{t) = (^v + 5X11 (t)) ^ + 6X^. (2.19) 
The corresponding linearized expression of the momentum reads: 

P = M, ^ = ~ (7 + ^h5X\^^ (v + 5X^ = + 5p, (2.20) 

yi-x -x 

where we introduced the zeroth-order Lorentz factor 7 = (1 — w ^)^^/^. The zeroth- 
order term is po = jMqV, and the longitudinal and transverse momentum fluctuations 
are given by: 

= 7Mg(l + i;V)<^^" = l^Mq6X\ 6pi = -fMg6X^. (2.21) 

It is convenient to separate the longitudinal and transverse components of the prop- 
agators, since as it will become clear in the next section, the off-diagonal components 
vanish: 

G^i(t) = Gll(t)^ + G^{t) (s^^ - (2.22) 



y2 \ y2 



and the corresponding decompositions for rj^^ and k*-' from (2.13). 

Inserting these expressions in equation (2.17), we find to zeroth order: 

^ = -vIpo, Po = iMv (2.23) 
and to first order in 5X the relativistic Langevin equations for position fluctuations: 

7X5^" = -ri\v)5X\\ + (e"(t)^"(^')) = n\\6{t-t'), (2.24) 

^M,5X^ = -rj^{v)5X^ + C^, {^{^^{1')) = K^5{t - t') (2.25) 

where the friction coefficients t]"'"*" are related to the coefficients 77^ by 



r]^ = 7Mgr/^, r/" = ^^M, ( 7]% + p^ 



p='yMqV 



(2.26) 



As we shall see in the following sections, the holographic prescription will directly 
compute the friction coefficients 77*-' and the diffusion coefficients k^^ appearing in 
equations (2.24-2.25). 
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Short-time solution: momentum broadening. For times shorter than the re- 
laxation time T£) we can treat the quark as travelling at a constant velocity 
V ( which is a good approximation in the case of a very heavy quark). In this 
regime^, one can write an approximate solution for equation (2.17), which describes 
a Brownian-like diffusion for momentum fluctuations. 

We start once again with equation (2.17), and linearize it (this time staying in 
momentum space) around a uniform trajectory p ~ Pqv/v + 6p. In the longitudinal 
and transverse directions we find the two equations: 



d6p^ 



d5p\\ 



Vd+P 



dp 



(2.27) 



PO 



where tj^^q = rf^ipo). The solution to these equations is straightforward: assuming 
initial conditions 5p{t = 0) = 0, it reads (notice that p-^ = 6p^): 



p^it) = / dt'e'^oA^'-'^^it'), 



Vd+P 



dt] 



D 



dp 



(2.28) 
(.2.29) 



PO 



From these solutions, we can compute the noise-average of the transverse and 
longitudinal momentum fluctuations 



((p^)2)= / dt' / dt" e'^oAi'+'"-^'\e{t')i^{t")) 



(2.30) 



{{p\^-p,f)=pl[l-e-^o,')\ I dt' I rft"e^"(*'+*"-2*)(e"(0e"(O)(2.31) 



Using the fact that {^{t')^(t")) = K6{t'—t"), and expanding to linear order in triD <^ 1, 
we arrive at the final result: 



{(p^f) = 2kH, ((Apll)^) = kH. 



(2.32) 



The first equation describes transverse momentum broadening, and it is typically 
parametrized in terms of the jet-quenching parameter q-^, 

q^ = ' ' = 2— 2.33 

^ vt V ^ ' 

The following sections, namely 4, 5, and 6 will be devoted to the calculation, in a 
5D holographic setup, of the Langevin correlators (2.4) and of the diffusion constants 
and At-*- appearing in equation (2.17). 



^As we are relying on the local form of the Langevin process, equation (2.17), we must still require 
time separations much larger than the auto-correlation time Tc. More explicitly, we consider time 
scales t such that Tc ^ t <^ tq- Therefore, consistency demands that td S> Tc- 
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2.2 The Langevin equation in the gravity dual picture 

As we have reviewed in the previous subsection, the memory kernel and the noise 
time-correlation function that govern the generalized Langevin equation (2.2) for an 
external quark, are given by appropriate real-time correlation functions of the force 
operator J^{t) over the ensemble that describes the medium. 

These correlation functions are precisely the kinds of objects one can compute 
in the gravity dual picture: one needs to identify the appropriate bulk field that 
couples to the boundary operator J^, then solve the bulk equations for this field with 
appropriate boundary conditions. 

As first discussed in [7, 8, 9], and as we will review in detail in Section 3, a probe 
heavy quark propagating through the plasma is described, in the gravity dual picture, 
by a probe string with an endpoint attached to a flavor brane, and extending into 
the bulk. The string endpoint moves along the quark trajectory and the rest of the 
string trails its endpoint extending in the holographic directions. The string world- 
sheet is described by the embedding coordinates X^{a, r) which, in the static gauge 
r = t, cr = r, reduce to the spacial components X(r, t), where r is the non-compact 
holographic direction. 

Using the trailing string picture for the heavy quark, the identification of the 
appropriate bulk field is straightforward: from equation (2.1) it is clear that the 
external source for the boundary field J^i{t) is nothing but the quark position X^{t), 
i.e. the boundary value for the string embedding X*(t,r). More precisely, for a 
heavy quark that follows an approximately uniform trajectory X^{t) = v^t + 6X^{t) 
the boundary coupling is of the form 

Scoupung = j dt6X\t)T,{t) , SX\t) = SX{n,t) (2.34) 

where 6X{rb, t) is the boundary value of the fluctuation in the trailing string around 
the classical profile. Therefore, the correlation functions (2.4) of the force opera- 
tors can be extracted, in the Gaussian approximation, by solving for the bulk linear 
fluctuations around the trailing string and using the appropriate holographic pre- 
scription. 

This calculation was first performed in [12, 14], for the AdS case. In this case, 
the world-sheet fluctuations propagate on a space-time with a metric of the form 
(3.30), with 

b{r) = i/r, f{r) = 1 - (7^Tr)^ (2.35) 

where T is the black-hole temperature. For a quark velocity v, the induced metric 
has a horizon at < < r^, with associated temperature Tg = T / The retarded 
correlator for the longitudinal and transverse components of J-" was determined using 
the prescription of [51]: 

r'Gi^AdS = GiAdS = -"^^O-dr^lsoundary (2-36) 
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where \l/(r, w) is a solution to the fluctuation equation with unit normahzation at 
the boundary and infaUing boundary conditions at the horizon, Q^^ = {27ri'^)~^H^'^ , 
with H^^ given by equation (3.35) speciahzed to the ArfS-Schwarzschild case. 

In order to compute both terms entering the Langevin equation, one needs also 
the symmetrized correlator, which gives the noise time-correlation function. In gen- 
eral, the relation between the retarded and symmetrized correlator depends on the 
statistical ensemble one is dealing with. For a black hole (as in the case of the 
induced world-sheet metric), the features of the ensemble can be obtained by con- 
necting the mode solutions along a Keldysh contour between the two boundaries of 
the maximally extended Kruskal diagram. This corresponds to obtaining a statistical 
ensemble by tracing over the degrees of freedom of one of the causally disconnected 
regions. In the context of AdS /CFT this idea was put forward in [13], which also 
provides a justification of the prescription (2.36) for the retarded propagator. 

We will not go into the details of this procedure, which can be found in [14, 15]. 
The crucial point is that the stationary statistical ensemble one obtains is a thermal 
ensemble at the temperature Tg = Tj Therefore, one can compute Gsym from 
the imaginary part of G/j, as in equation (2.14), with the substitution T — )■ Tj^. 

Notice that the retarded Green's functions compute, through equations (2.13), 
the coefficients appearing in the Langevin equations for the fiuctuations equa- 
tion (2.24-2.25). In particular, to extract the coefficient "(}\) must divide G^\lS) 
by an extra factor 7^ with respect to the corresponding result for 77^. 

Computing the zero-frequency limit of the retarded correlators, the resulting 
Langevin diffusion coefficients are found to be [14, 15]: 

7-2 J =^^ = J^^^T\ (2.37) 
and the friction coefficients reproduce the classical drag force calculation [9], 

The diffusion and friction coefficients indeed satisfy an Einstein relation appropriate 
for the temperature T^: 

n^/ri^ = /r/ll = 2MT„ (2.39) 

where 77"*" and 77" are related to rjn by equation (2.26). Notice that, in the conformal 
case, 77D is momentum-independent. 

The approach taken in [14] derives the Langevin propagators by using the stan- 
dard holographic prescription for real-time correlators. An alternative procedure, 
giving the same result, was adopted in [16, 17] for the non-relativistic case, and later 
in [18] for the general relativistic case. These authors performed a direct derivation 
of the Langevin equation: starting from the trailing string fluctuations in the bulk. 
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and integrating them out, they showed exphcitly that one arrives at equations hke 
(2.24-2.25) for the boundary fluctuations, with coefficients given by the formulae 
previously found in [12, 14]. Furthermore, they showed that the same result can be 
obtained by integrating out only a strip between the world-sheet horizon and the 
stretched horizon -|- e: this gives a picture of the stochastic behavior of the string 
fluctuations as originating from the world-sheet horizon. In this work we will not 
follow explicitly this road, but rather rely on the holographic computation of the 
Langevin correlators. Nevertheless, by using the general formalism developed [45] 
we will compute directly the transport coefficients. 

3. 5D non-conformal backgrounds for Langevin holography 

In this section, we present the background material for the holographic computation 
of the Langevin correlators that we carry out in the Section 4. As a bulk geom- 
etry we take a general, five-dimensional, asymptotically AdS black-hole, dual to a 
non-conformal deconfined plasma. These geometries arise generically as solutions in 
appropriate Einstein-dilaton theories in five- dimensions [35, 36]. 

A heavy external quark moving through the plasma at temperature T can be 
described by a string whose endpoint at the boundary follows the quark's trajectory 
[7]-[31]. The string extends into the bulk, whose geometry is the dual black hole 
background with appropriate temperature T. Once the motion of the endpoint at the 
boundary is specified, one can find the trailing string solution through the geodesic 
equation: the momentum flow along the string is dual to the drag force experienced 
by the quark moving through the plasma. 

The fluctuations of the string world-sheet around the geodesic solution are holo- 
graphically dual to the stochastic forces felt by the quark due to its interaction with 
the medium. Their effect to leading (Gaussian) order is that of stochastic noise act- 
ing on the quark, resulting in a Langevin-type diffusion with its associated transport 
coefficients (diffusion constants). 

3.1 5D Einstein-Dilaton black holes 

We shall consider the dynamics of a probe string in a general 5D black hole geometry, 
with a string frame metric: 



f{r)dt^ + dx^dxi 



(3.1) 



We assume there is an asymptotically AdS region r — t- where 



log h{r) ~ — log - + subleading, /(r) ~ 1 -|- 0{r'^) 



r 



(3.2) 



r 
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and a horizon a.t r = where f{rh) = 0, and f'{rh) and b{rh) remain finite. The 
black hole temperature is given by: 

4vrT=-/'(r,). (3.3) 

We make no particular assumptions on the subleading terms in equation (3.2). 
In case these subleading terms actually vanish as r — )■ 0, then the metric is AdS in 
the usual sense. 

The black holes of the type (3.1) arise, in particular, as solutions of a large class 
of 5-dimensional Einstein-dilaton models, described by the Einstein-frame metric gj^^ 
and a scalar field A with the action: 

(3.4) 



S = -MlNl j 



In the holographic interpretation of these models, the scalar A is dual to the running 
coupling \t of the four-dimensional gauge theory. This is the class of models we will 
have in mind, although the results of this work apply to any 5D theory that admits 
solutions such as (3.1). 

For an appropriate choice of the potential ^(A), the models with action (3.4) 
provide a good holographic dual to large- A'^c 4-dimensional pure Yang-Mills theory, at 
zero and finite temperature [32]- [38]. The potential should have a regular expansion 
as A — 0, with 

F(A)~^(l + i;oA + ...) . (3.5) 

Furthermore, linear confinement in the IR requires that, at large A, V{X) grows at 
least as fast as X^^^. With these requirements, 

1. The solutions in the Einstein frame are an asymptotically AdS metric, with 
AdS length and a non-trivial profile A(r); 

2. There is a first order Hawking-page phase transition with a non-zero critical 
temperature Tc. 

For a short review of the main features of these models, the reader is referred to [6]. 

We will be interested in the fiuctuations of a probe string in the 5D black hole 
geometry for T > (corresponding to the deconfined phase). The black hole solu- 
tions in the string frame have the form (3.1), with string frame metric gfj,u = \'^^^g^y 
In the UV region r = 0, we therefore have: 

h{r) ~ -\''\r\ f{r) ~ 1 - § r^ A(r) ~ (3.6) 
r i-^ 8f log r 

where the constant C depends on the thermodynamic quantities that characterize 
the black hole, and it can be expressed in terms of the temperature T and entropy 
density s, [36]: 

C 457r2 _ ^ ^ 
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The expressions (3.6) are corrected by terms of 0{X), which are neghgible near 
the boundary. 



3.2 Classical trailing string and the drag force 

Before going into the details of the world-sheet fluctuations, we review the calculation 
of the unperturbed trailing solution, that was discussed in [7, 9] for pure AdS black 
holes, and generalized in [11] for black holes in 5D Einstein-Dilaton theories. In this 
subsection we review the setup and the results of [11]. 

We consider an (external) heavy quark moving through an infinite volume of 
gluon plasma with a fixed velocity v a.t a. finite temperature T. In the dual picture, 
this is described by a classical "trailing" string with an endpoint on the UV boundary 
moving at constant velocity v. 

The world-sheet of the string is described by the Nambu-Goto action'' , 

Sng = j c^Vv^-det , gap = g^.udaX^dpX" , | ^ ^'^ 

' (3.8) 

where gfj,^ are the components of the bulk metric in the string frame^, 



ds' = h\r) 



- f{r)dr + dx'dx' 



1,2,3. (3.9) 



.fir) 

The ansatz for the classical trailing string is [9] , 

X^ = vt + ^{z), X^ = X^ = 0, (3.10) 
and along with the gauge choice 

e = t, e = r (3.11) 

leads to the induced metric: 

s»,=''^(o('':,5'-' ,5,<'-.',.V (3.12) 



and the corresponding action: 



S = ~ I dtdr h\r)Jl - + /(r)r^(r) . (3.13) 



^Throughout the paper, we will denote 5D coordinates by ji^v . . world-sheet coordinates by 
a,(3..., and boundary spatial coordinates by i,j . . .. Indices i,j ... in the boundary theory are 
raised and lowered with metric Sij , so we will make no distinctions between upper and lower indices 
as far as boundary tensors are concerned. 

^Unless otherwise stated, b{r) will always denote the scale factor in the string frame. This is a 
shght change of notation with respect to our previous papers [32, 33, 35, 36, 38], where the same 
quantity was denoted bs{r), but it is justified since most expressions in this work are simpler in the 
string frame. 
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Since S does not depend on ^ but only its derivative, the conjugate momentum tt^ 
is conserved, 



where the final expression is obtained by evaluating it at the point r = r^, defined 
by 

f{rs)=v^. (3.15) 

For an infinitely massive quark, the string endpoint is the boundary, r = 0. 
For a quark of finite mass Mq, the endpoint should be located at a position rq in 
the interior, as discussed in detail in [11]. This puts an upper bound on the quark 
velocity v, since the trailing string picture fails when < tq. At this point, the 
fiavor brane dynamics should become important. 

The Drag Force The drag force on the quark can be determined by calculating 
the momentum that is lost by fiowing from the string to the horizon, which results 
in: 

i^drag = TT^ = (3.16) 

where we have replaced /(r^) by v"^ in the last equality. 

One defines the momentum friction coefficient r/^ as the characteristic attenua- 
tion constant for the momentum of a quark of mass Mq-. 

i^drag = ^ = -VdP, , P = MgV-f (3.17) 

where 7 = (1 — t;^)"^/^ is the relativistic contraction factor. With this definition we 
obtain: 

In the conformal case, rjD is independent of p, 

conf TtVXT^ 

where A = {i/ig)'^ is the fixed 't Hooft coupling of = 4 sYM. This is not anymore 
so in the general case, where r]D is momentum dependent. 

The world-sheet black hole. The coordinate value r = is a horizon for the 
induced world-sheet metric. In order to ascertain this, we can invert equation (3.14) 
to obtain ^'(r) in the form: 
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We may now change coordinates to diagonalize the induced metric, by means of the 
reparametrization: 

/ - ^ /(r) V(/(r) -t;2)(64/_C2) 



In these coordinates the induced metric is 
ds^ = b 



(3.22) 



The coefficient of c/r^ vanishes at , so this point corresponds to a world-sheet black- 
hole horizon. Since /(r) runs between and 1 as < r < r/j, by definition (3.15) the 
world-sheet horizon is always outside the bulk black hole horizon, and it coincides 
with it only in the limit — )■ 0. In the opposite limit, f — )■ 1, asymptotes to the 
boundary r = 0. 

The Hawking temperature associated to the black hole metric (3.22) is found 
as usual, by expanding around r = and demanding regularity of the Euclidean 
geometry. The resulting temperature is: 




/(r.)f(r. 



4y(r.) ^ firs) 
b{rs) firs) 



(3.23) 



In the conformal limit, where the dilaton is constant and the background solution 
reduces to AdS'-Schwarzschild, the world-sheet temperature and horizon position are 
simply given by: 

= ^, r^f = (3.24) 

More generally, in the ultra-relativistic limit ^ 1, one can express in terms 
of thermodynamic quantities. In this limit approaches the boundary r = 0, and in 
this region the geometry approaches that of AtiS'-Schwarzschild, equations (3.6,3.7). 
Therefore, from the definition (3.15) we obtain: 

3.3 Fluctuations of the trailing string 

We now proceed to study the quadratic fluctuations around the classical trailing 
string solution reviewed in the previous section. This analysis was performed in the 
AdS black hole background in [14]. Here, we extend it to the general 5D background 
(3.1). 

We continue to work in the static gauge = t, C,^ = r, but we will allow for a 
more general ansatz for the embedding coordinates: 

=vt + ^{r) + 6X\^{r,t) , = 6X\r,t) , X^ = SX^{r,t). (3.26) 
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We will treat the quantities 5X* as perturbations around the background solution 
(3.10). 

The Nambu-Goto action (3.8) is now given by: 



S 



NG 



2tiP 



where 



f 



X''X'' 1 , grt = b'^ X' X' 



(3.27) 



(3.28) 



where a dot and a prime represent derivatives w.r.t. t and r respectively. 

Expanding the Nambu-Goto action in 6X^ around the classical solution (3.20) 
we obtain, to quadratic order: 



^2 



27r£2 



dtdr G"'^ 



^dJX^^d^dX^^ + ^Y1 9JX'dp5X' 



i=2 



where 



G 



al3 



and we have defined: 



Z3 



Z = h^ 



b^f-C^' 



(3.29) 



(3.30) 



(3.31) 



Note that det{G°'^) = —b^/Z^ and that in the J\f = A case Z = \/l — v"^ is a constant. 
In terms of the induced world-sheet metric (3.12), we obtain 



= Z-^fc^^"'^, V- det g = b^Z . 
We may therefore rewrite the action as 



S, 



._L_ f dtdr^l ^ - det g g'^^ ^dJX\\d^6X^\ + Y,dJX'dis6X 



(3.32) 



(3.33) 



To simplify the action, we change coordinates to diagonalize the induced metric, as in 
the previous subsection. By a reparametrization of the world-sheet time coordinate 
as in (3.21), the new induced metric is (3.22), and the action read: 



^2 



27r£2 



drdr -H'^^ 
2 



1 ^ 

-dJxWdpdX'^^ +^dJX'dp5X' 



Z2 



with 



i=2 







^/{f-v^Wf-c^) 



(3.34) 



(3.35) 
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Equations (3.29) and (3.34) show that the longitudinal fluctuations (i.e. those parallel 
to the direction of the unperturbed trailing string motion), namely 5X\ and the 
transverse fluctuation 6X'^ and 6X^, have different kinetic terms, as the effective 
two-dimensional metrics they are sensitive to differ by a factor Z^. From now on, 
we will denote the longitudinal fluctuations as SX^^ and the transverse fluctuations 
as 6X^. 

From equation (3.22) one can immediately derive the field equations satisfied by 
the fluctuations: 

da Z-^ H^'^d^SX^^ = , 9„ H^'^dpdX^ = 0. (3.36) 

For a harmonic ansatz of the form (5X*(r, r) = e^'^'^6X^{r,u), equations (3.36) 
become: 



2 hi 



dr Jif -v^Wf -C^ driSX^) + , 5X^ = (3.37) 

dr l^Vif -^'^Wf -C^) dr (SX^] + =^^^^=6X^\ = (3.38) 

In the next sections we will compute the Langevin correlation functions from these 
fluctuation equations and extract the diffusion constants and the spectral densities 
from them. 

We note however that the diffusion constants can also be read-off directly from 
the quadratic action (3.34) by using the method of the membrane paradigm as ex- 
plained in section 4.2. 



4. Holographic computation of Langevin correlators 
4.1 The Green's functions 

From the discussion in the previous section, it emerges that in a 4D theory with a 
5D gravity dual we can compute the Langevin correlators holographically, from the 
classical solutions for the fluctuations of the trailing string. As we have observed 
in section 3, these fluctuations behave as free fields propagating on a 2D black-hole 
background, whose metric is essentially the induced metric on the bulk trailing string, 
equation (3.22). The asymptotic form and the causal structure of this black hole are 
exactly the same as the one for the trailing string embedded in an AdS black hole. 
As a consequence, the results of [14, 15, 18] discussed in the previous sections imme- 
diately generalize to the more general metric (3.1): following the Keldysh contour 
in the extended Kruskal diagram of the black hole, one finds a thermal spectrum 
of transverse and longitudinal fluctuations with effective temperature Tg, given in 
equation (3.23). 
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In this ensemble, the symmetrized and retarded Green's functions obey the re- 
lation (2.14), with T = Tg. Therefore, one can obtain both the memory kernel and 
the noise correlator entering equations (2.2-2.3) from the knowledge of the retarded 
Green's function. 

From the structure of the action for the fluctuations, equation (3.33), one can 
observe that there are essentially two types of retarded correlators, G|j(aj) and Gj^lu), 
for the longitudinal and transverse fluctuations. We introduce a notation similar to 
[14] and define 

1 1 fJClP 

pap _ J- TTo/3 paP _ -L -n IA^\ 

= ' ^11 = 2^^' ^^-^^ 

where H""^ is defined in equation (3.35). 

The holographic prescription for the retarded correlator, computed with the 
diagonal induced metric (3.22) is given by: 

Gnioj) = - [^Ur, Uj)g^''dr^n{r, U;)]boundary • (4-2) 

Here ^R{r,u) denotes collectively the fluctuations 6X\ 6X-^, solutions of equations 
(3.37-3.38) with the appropriate boundary conditions, i.e. unit normalization at the 
boundary, and infalling conditions at the world-sheet horizon (as we discuss more 
extensively below) and the factor Q"' is the appropriate one from equation (4.1) 

The expression in equation (4.2) must be evaluated at the boundary of the trail- 
ing string world-sheet. In the case of an infinitely massive quark, the string is at- 
tached at the AdS boundary at r = (when needed, in order to keep quantities 
finite, we introduce a cut-off boundary at r = e). In case we want to keep the quark 
mass finite, the trailing string is attached to a point rg, which is determined by 
demanding that the free energy of a static string in the T = background gives the 
mass of the quark: 

M, = ^ £ h{rfdr (4.3) 

(here r^, is the point at which h{r) reaches its minimum, see the discussion in [11]). 

Next, we discuss in greater detail the boundary conditions for the rs. The 
solutions to the fluctuation equations (3.37)-(3.38) share the same asymptotics both 
for the transverse and longitudinal components (since Z{r) asymptotes to a constant 
both at the horizon and at the boundary, where the equations have singular points). 
At the world-sheet horizon r — equations (3.37)-(3.38) both take the form 

dl^ + ^^M+(-^ 1") 5^ = 0, (4.4) 

\r-rs\ \'i'KTs\r - r s\ J 

so that the solutions near the horizon behave as 

^{r,uj) ^ {r-s-r)^*^' ^ (4.5) 
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The + sign in the exponent corresponds to a wave which is outgoing with respect to 
the world-sheet horizon, while the — sign characterizes an in-falling wave. 

Near the boundary r — )■ both transverse and longitudinal fluctuations have to 
solve the following equation: 

d^<i/{r,u) -(-- ~) dr^{r,u)+^W<i/{r,u) = , (4.6) 
\r 6 A J 

As long as r A'/ A ^ 1 ^ the two independent solutions are a normalizable mode and 
a non-normalizable mode, 

^ ~ + Cy\-^/^. (4.7) 

According to the standard prescription [51], the appropriate boundary conditions 
for the wave functions in the expression (4.2) for the retarded correlator are the 
in-falling behavior at the world-sheet horizon with the condition ^ R{r) = 1 at the 
boundary: 

^nir^uj) ^^!h {rs-r)~^s r r^. (4.9) 

where is a constant. 

Given the wave-function, obeying the near-boundary and near-horizon asymp- 
totics specified by equations (4.8)-(4.9), we can extract the propagator from equation. 
(4.2). Below, we separately discuss the features of the real and imaginary parts of 
the retarded Green's functions, and the associated spectral densities. 

Real part of the retarded correlators. The real part of the correlator (4.2) 
suffers from ambiguities related to the possibility of adding boundary counterterms 
to the action (3.8). This was discussed e.g. in [52] in the context of the calculation 
of 4-dimensional transport coefficients. 

The ambiguities in the propagator usual, associated to UV-divergences in 

the on-shell action, that arise when we try to evaluate it in the full AdS space-time. 
To obtain a finite result, we must consider the action on a regularized space-time, 
with boundary at r = e rather than r = 0. Then, once the divergences in the limit 
e — >■ are identified, one can add counterterms to subtract them and obtain a finite 
limit as e — 0. 

The essence of holographic renormalization is that these counterterms are local, 
covariant boundary terms. As we show in appendix A, we only need a single boundary 

^The condition rA'/A ^ 1 is realized in particular in the case of logarithmic running: in that 
case, for small r, rA'/A ^ A << 1. In the case where A is dual to a relevant operator it is also valid 
a fortiori since rA'/A ~ Ar, with A = mm(A, 4 — A). 
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counterterm to regularize the action, and this is given by the boundary-covariant 
point-particle action: 

S^ount = AM(e) J dt ^X>^X^ (4.10) 

The same UV divergences appear in the real part of the propagator, if we try to 
compute it naively from equation (4.2). In fact, the expression (4.2) is nothing but 
the unrenormalized on-shell action, as can be easily observed by integrating equation 
(3.34) by parts and using the field equations (3.36). Therefore, as a consequence of 
the analysis of the on-shell action carried out in Appendix A the divergent parts of 
the transverse and longitudinal Green's functions are: 

«-«)'^'^'^^4-^ ,,,, 

This result can also be explicitly derived from the explicit form of the wave- 
functions, close to the boundary. As we will show in Subsection 4.4, (see eq. (??), 
and for more a more detailed derivation. Appendix C) the solution of eq. (4.6) is, 

^uv{i^) = [cos(7Ci;r) + {'jur) sin(7Ci;r)] + CyX~^^^{r) [{'jur) cos{'-fur) — sin(7c<;r)] 

(4.12) 

This solution generalizes eq. (4.7) for any finite u, and it is valid in the near- 
boundary region, i.e. for 'yur <^ 1, and A(r) ^ 1. The value of the coefficient Cs 
of the leading term is fixed to ensure unit normalization at r = e (we are keeping in 
mind that we will take the e — )■ limit at the end). 

Evaluating the real part of (4.2) at r = e with the wave-function given by (4.12), 
we find the following divergent term for the transverse and longitudinal components: 



-2„„.,,..A'"(«).,... _ 



Re ~ 7-2Re ~ —-fu' 



1 + le^co' + 0(6^) 



e^O (4.13) 



The divergence is purely in the u"^ term, and the coefficient agrees with the result we 
found from the on-shell action (4.11). Notice that the second term in equation (4.12), 
proportional to C^,, starts at 0{e^u^), so it does not contribute to the divergent part 
of the propagator. 

To eliminate the divergence, and obtain a finite result, we must add the contri- 
bution from the boundary counterterm. However, different results can arise due to 
different choices for the finite contributions included in the counterterm, i.e. different 
subtraction schemes. As discussed in appendix A, in our case this ambiguity reduces 
to a term of the form 6Gr{uj) = 6m u"^, which can be reabsorbed in the renormaliza- 
tion of the quark mass. We are going to use a minimal scheme, and fix the coefficient 
of the counterterm action (4.10) to be: 
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This choice exactly subtracts the divergences (4.11) (see Appendix A), and moreover 
removes all (finite or infinite) 0(a;^) terms in the large-cu behavior of the propagator. 

Once we subtract the divergence in the minimal scheme, and we take the limit 
e — )■ 0, the right hand side of equation (4.13) vanishes. This means that, in this 
scheme. Re Gfi{uj) — )■ as a; — )■ oo. Moreover, as the only ambiguity in Re Gr is 
proportional to u"^ (see Appendix A), we conclude that in any other scheme the real 
part of the Langevin Green's function grows as for large co. 

As a final remark, the previous discussion only applies if we consider the quark 
mass to be infinite. In the case of a finite mass, the trailing string is attached at a 
radial point rg > 0, the cut-off is physical, and the result is not divergent. However, 
one should still specify the finite boundary term included in the action in order to 
arrive at an unambiguous result. 

Imaginary part, and the symmetric correlators. Unlike the real part, the 
imaginary part of the retarded correlator does not suffer from ambiguities. One of 
the reasons is that it is proportional to a conserved quantity, which can be shown to 
be finite at the horizon. In fact, we can write Im Gji in the form: 

Im Gr{u) = -ie;'^^v|/^S^v|/^ = -r . (4.15) 

Here is a conserved current — this follows directly from the equations for the fiuc- 
tuations, equations (3.36) — hence the imaginary part of the retarded correlator can 
be analytically evaluated at any r, not necessarily at the boundary. It is convenient 
to evaluate it at the horizon. From the definitions (4.1) and (3.35), we find that, in 
the near-horizon limit : 

QY ^ Z\r,)gif ~ (47rr,)6(r,)2(r, - r), r ^ r,. (4.16) 

Inserting this expression in equation (4.15), and using the expression (4.9) for \E'(r), 
we find: 

imG^(c.) = -^|v^^(c.)rc., i-4M = -^Jpl|;ryl^llMr^> (4-i7) 

where \E'/i is the coefficient of the in-falling wave-function, see equation (4.9). 

From the imaginary part of Gr{uj) we can immediately extract the symmetrized 
correlator Gsym{^), i-e. the generalized Langevin noise time-correlation function: as 
discussed at the beginning of this section, due to the thermal nature of the world- 
sheet fiuctuations, Gsymi^jj) is related to Im Gr{(jj) by the following equation, which 
generalizes the analogous equations in [14, 15, 18]: 

Gsym{i^) = -coth Im Gr{u) . (4.18) 
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The spectral densities. The spectral densities associated with the Langevin dy- 
namics are defined by 

p«H = --Im G^oo), (4.19) 

TT 

PtyM = -lo^^mi^) = -COth (^^^ p^iu), (4.20) 

where a =_L, || , and we have used equation (4.18) in the second line. The imaginary 
part of Gfi is given by the fiux for the perpendicular and the parallel components 
by the formula (4.15). We will give an analytical estimation of the large-frequency 
behavior of the spectral densities in Subsection 4.4, using a WKB method; in Section 
6 we will use numerical methods to obtain the full functional dependence of p{u) in 
a concrete model. 



4.2 The membrane paradigm 

Here we introduce an alternative method to calculate both the diffusion constants 
and the spectral densities that goes under the name of the membrane paradigm [45] . 
We apply this method to obtain the spectral densities in the next section and to 
obtain the diffusion constants, in section 5.3. 

In [45], it was established that the transport coefficients associated with generic 
massless fluctuations can be read off directly from their effective coupling in the 
action, evaluated at the horizon. For an arbitrary massless fluctuation (p with an 
action 

S2 = -IJ d'^xdr^J- detgQ{r) (4.21) 
the transport coefficient associated with the retarded Green's function is given by. 



\mGR{uj,k) -detg 
Xr = ~ hm = lim Q[r.^ ' 



(4.22) 



where Q is the effective coupling of the fluctuation deflned in (4.21). We refer the 
reader to [45] for a derivation of this formula, noting that here we apply the same 
idea to the world-sheet black-hole rather than the bulk black-hole as in [45]. 
In our case, comparison of (4.21) with (3.33) yields. 

Once we know these effective couplings, we can immediately write down the diffusion 
constants. Therefore, the method provides a very efficient and fast way of computing 
the latter. This is done in section 5.3. 
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4.2.1 A differential equation for spectral densities 

Although the method of the membrane paradigm is most effective in the low fre- 
quency limit, where one can read off transport coefficients directly from the Q's, it 
is still a convenient method for arbitrary uj where one has to use the flow equations 
[45]. We define the canonical momentum associated with in (4.21) with respect to 
foliations in r, as 

n{r,u,k) = Q{r) ^- det ^/"9,(/>(r, k), (4.24) 

where we performed the Fourier transform on the 4D space-time. We also define the 
"r-dependent" response function 

X = . "^^'"^'^j . (4.25) 
a;0(r, A;) 

Using the general AdS/CFT relation for the retarded Green function, 

G«(c.,fc) = lim^|^?^, (4.26) 

''^0 0(r,a;, k) 



we find that, 



p(w, k) = --Im Gr{uj, k) = - limRex- (4.27) 

7r vr r--s>o 



From the equations of motion, one derives a first-order equation for x, [45 

i 9rr ( 1 -2 



^,X = ^u:,^—[7^,r-Q{r)]. (4.28) 
9tt \Q{r) 



Here, the effective coupling Q is given by (4.23) and, 

Z(r) 




fir) - 



(4.29) 



For regularity at r = one should require that X Qi^s) as r — r^. One solves 
(4.28) with this boundary condition at the horizon, and determines x on the boundary 
r = Tfj. The spectral density associated with the symmetric Green's function is given 
by, 

Psymiuj) = ^ COth (^^^ Rf^Xijb, w). (4.30) 

As we show in section 4.4, psym is divergent as w — )■ oo. This is similar to the familiar 
short-distance divergences of correlators of quantum field theory. 



10 



In the rest of this section, we set fc = 0. 
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4.3 Universal results for the spectral densities 

The membrane paradigm allows us to obtain interesting relations concerning the 
spectral densities in certain limits. We first note that, by employing the equation 
of motion for that follows from (4.21) and (4.28) one can show that |0pRe x is 
independent of r: 

|</.(r)|2Rex(r) = iV(a;). (4.31) 

Through equation (4.27), this is of course equivalent to the fact that the imaginary 
part of the Green function is proportional to the conserved flux, hence it is constant. 

We may evaluate (4.31) in the two limits r = rt and r = r^. We flrst evaluate it 
at the horizon: the in-falling condition is, 

0(r,w) -> C^(r„w)(l -r/r,)"3^, r ^ r,. (4.32) 

Therefore, |0(rs)p = |C/jp. On the other hand, by the boundary condition of (4.28), 
xi^s) = Q{i^s)- Therefore we obtain, N{u) = Ch{(jJ,rs)Q{rs). 

Secondly, we evaluate (4.31) at the boundary ri, using (4.27) and we obtain 
A^(c<;) = 7ip{u)/u. Hence we have, 

p{u) = -Q{rs)\Ch\'{rs,u). (4.33) 
vr 

Now, we consider the special limit — > r^, while keeping u flnite, or more precisely 
we consider uTs <^ 1. If we think of u as being flxed, this limit can be attained in 
two ways: 

• either by sending T — )■ oo, so that the black hole horizon, and consequently 
the world-sheet horizon, are pushed to the AdS boundary; 

• or by keeping T flnite, and sending v —> 1. In this case the black hole temper- 
ature is flxed, but the string is ultra-relativistic. 

Therefore, the regime uvg ^ 1 corresponds to a UV limit for the background quan- 
tities T and v. 

On the other hand, for a given flxed u, as approaches the boundary Ch{oj) 
becomes independent of u and approaches unity. This is because of the unit nor- 
malization of at the boundary, 0(rb) = 1. Therefore, (4.33) simplifles and one 
obtains, 

p{uj) ^ -g(r,), ooTs < 1 (4.34) 

TT 

We can make this expression more explicit by taking the near-boundary approximate 
expressions, valid for r<, ~ 0, in (4.23): we approximate the metric functions as in 
equation (3.6), and consequently, from equation (3.31), Z{rs) ~ 7"^. These approxi- 
mations are valid up to terms of 0{X), which signal the departure from conformality 
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of the UV geometry. Using these approximations we find: 

Q±irs) ^ 7^'Q||W ^ ^^%^[l + 0(A(r,))l. (4.35) 

Finally, can also be expressed purely in terms of the boundary quantities in 
the same limit, using the near-boundary approximation (3.25). 

The conclusion of the previous discussion is that, in the high energy limit one 
obtains universal expressions for the spectral densities, where they become linear in 
frequency. 

This universal behavior is to be expected based on the consideration that, for 
fixed u, both limits f — ?■ 1 and T — oo correspond to the limit — )■ 0. In this 
regime the equation governing the string fiuctuation become essentially the equation 
one finds in pure AdS, close to the boundary. This case is discussed explicitly in 
Appendix D, from which one concludes (see equation (D.l)) that the wave-functions 
depend on u and Tg only through the combination u/Tg oc uvg. Thus, for fixed u 
and small r<,, i.e. for uvg <C 1, the spectral density is approximated by the linear 
term in the expansion in uVg. 

4.4 The WKB approximation at large frequency 

The WKB approximation can be used to obtain the high-frequency behavior of the 
Langevin spectral densities. Here we summarize the method and present the results. 
The detailed derivation is given in appendix C. 

By a rescaling of the wave-function, the fiuctuation equations (3.37)-(3.38) can 
be put in a Schrodinger-like form, and the large u solution can be obtained by an 
adaptation of the WKB method. This method has been applied to the case of shear 
perturbations in AdS-Schwarzschild black-hole in [46]. 

The fiuctuation equations (3.37)-(3.38) can be put in the Schrodinger form 

-r + K(r)^ = 0, Vsir) = + ^(log7^)" + ^(log7^)'^ (4.36) 

defining the wave function ^ = 71^"^ and R = •\/(/ — v^){b'^f — C^), where 

For large frequency (compared to r~^), we can approximate the potential over 
essentially the entire range of r by the expression: 

Vs{r) ~ > 1, rtp<^r< r„ (4.38) 

where r^p denotes the classical turning point, V{rtp) = 0. The range of r for which 
equation (4.38) is valid is the classically allowed region. For large uvg, the turning 
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point is approximately r^p ~ v^/(7i^) ^ ^g. Therefore, the classically allowed 
region covers almost all the range < r < r^, but for a small region close to the 
boundary (that includes the turning point) where the approximation (4.38) breaks 
down. Finally, and most importantly, for uvs ^ 1 there is always an overlap region 
in which the boundary asymptotics of the Schrodinger equation and equation (4.38) 
are both valid , since r^p ^ r^. This will allow the matching of the WKB solution 
and the boundary solution. 

In order to obtain a large-w analytic approximation for p{uj), it is necessary to 
specify whether we consider the quark mass as infinite, or we are working with a 
finite but large mass Mq. 

Infinitely massive quarks. In this case the endpoint of the string is attached to 
the AdS boundary r = 0, and this is where we should set the normalization of the 
wave-functions. 

The WKB computation proceeds along the steps detailed in appendix C. The 
WKB solution is matched with the horizon and boundary asymptotics, for the wave 
function, determining the coefficients appearing in equation (4.17) 

, j . (4.39) 

b{rs) [lZ{rs), II 

Inserting these expressions in equation (4.17) and using the definition (4.19), we 
determine p{oj), in the limit u ^ I/t^: 

P±{uj) ^ 7 Vll(c^) - K(^^ (4-40) 

Here Xtp = \{rtp) where rtp ^ ^/2/(u'j) is the classical turning point, as discussed 
in Appendix C. For very large u, the dilaton profile is can be approximated as in 
equation (3.6), 



Xtp ^ bQ ^ log 



(4.41) 



Finite mass quarks. The computation for the finite mass quarks follows the same 
steps as for the infinitely massive case, except that the boundary normalization of 
the wave function \l/(r) = 1 has to be imposed at the cutoff r = tq (determined by 
equation (4.3)), rather than at the proper boundary r = 0. As Mq becomes large 
(with respect to the UV scale A), rg ~ ^/Mq. 

The presence of a finite cutoff at tq implies some subtleties in the matching of 
the WKB solution with the boundary solution, as explained in Appendix C. 

To give an explicit result, we must distinguish two regimes which, for fixed u/Mg, 
correspond to small and large velocities. 
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Small velocities. If tq <^ I/7W the finite mass spectral densities at large 
frequencies behave like the infinite mass ones, except for {'^ujrqY corrections. 
More explicitly, we obtain 



3 1)2 



TT 



4 

2 f2 Ki"^) 



1 + (jurQf + 



7^Q 



^« dr' 



n{r') 



(sin(7c<;rQ) - 'yurq cosi'^urq))" 



(4.42) 



These expressions are valid in the regime where uvg ^ 1, but at the same 
time 'yoJVQ < 1, i.e. for small velocities, given a fixed (large) quark mass and 
a given frequency. On the other hand, for fixed frequency and velocity but for 
Mq — > 00, vq — 7- 0, the r.h.s. asymptotes to unity, and we recover the infinite 
quark mass expressions (4.40). 



Large velocities. Analogously, in the limit where 
spectral functions read 



> Tip ~ A/2/7CU, the 



7 P\\ 



1 + {-furQY 



(4.43) 



Hence, the difference with respect to the infinite mass result in this case is 
that A^p^ — which is cu-dependent — is replaced by Ag'^ — which, instead, is 
w-independent. 



We note that the large uj behavior for finite mass, both in the large and in the 
small cutoff regimes, changes with respect to the infinite mass case and becomes linear 
rather than cubic, due to the extra {'^ujrQY term in equations (4.42)-(4.43). This 
extra term comes from the fact that the solution has a subleading linear dependence 
on r, which is negligible in the infinite mass case, but enters the expression of the 
spectral function, p ~ v]>*v]/'^ in the finite mass case, giving it a quadratic dependence 

on Tq. 



5. Langevin diffusion constants 

The correlators and spectral densities discussed in the previous section are the build- 
ing blocks of the generalized Langevin equation, (2.2). Now, we will focus on the 
long-time limit, discussed in Section 2, in which equation (2.2) reduces to the local 
form (2.12), in which only the cj-independent friction and diffusion coefficients, 77 and 
K, appear. They are given in equation (2.13) in terms of the zero-frequency limit of 
the Langevin Green's functions. 
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Therefore, in this Section we consider the zero-frequency hmit of the Green's 
functions constructed in Section 4. This will allow us to give the analytic results for 
the diffusion constants, both from the direct evaluation of the correlators, and using 
the membrane paradigm. 

5.1 Diffusion constants via the retarded correlator 

The diffusion constant is defined in terms of the symmetric correlator Gsym (see 
equation (2.13): 

K = lim Gsvm = — 2Ts lim ^ , (5.1) 

where in the second equality we have used the a; — j- limit of equation (4.18). 

The small frequency limit of the symmetric correlator can be evaluated ana- 
lytically since we can determine the boundary-to-bulk wave function in this limit 
and discard higher orders in the evaluation of (5.1). More precisely, we write the 
small frequency limit of the horizon asymptotics of the \&r's. Given the in-falling 
boundary condition (4.9) we obtain 

^^i(r, u) = ^h{rs - r)"3^ ~ 

This solution can be connected to the boundary asymptotics by the exact solution 
of the fluctuation equations (3.37)-(3.38) at u = which reads "^r^t, 0) = 1 once we 
impose the appropriate boundary condition \E'/j(rb) = 1 (see Appendix B for details). 
On the other hand, equation (5.2) reduces to = in the strict a; = limit. 
Therefore, the near-horizon solution at small frequencies is entirely determined by 
equation (5.2) and the match to the boundary solution which yields = 1 (both 
for transverse and longitudinal modes). 

Furthermore, expanding for a; ^ 1, we may also show (see Appendix B) that 
the solution for all values of r and small frequencies is given by 

^R = ir,-r)-^ [l + 0{u)] . (5.3) 

Now we may substitute the solution (5.3) in the expression for k, using (5.1) and 
evaluating the current at the horizon. 

Infinitely Massive Quarks. For infinite mass quarks, the boundary is located 
at r = and the appropriate boundary-to-bulk wave function is given by 
equation(5.3) at small frequencies. 

To compute the diffusion constants (5.1), we evaluate J^, defined in equation 
(4.15), at the radius value r = — e and then let e — (since J'' is conserved it can 
be evaluated at any radius and not necessarily at the boundary, where subleading 
0{uj) terms in (5.3) would contribute). This allows to neglect the sub-leading terms 



4^ 



log |r 



(5.2) 
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in the solution (5.3). For the longitudinal component we also need the near-horizon 
limit of Z{r), which can be easily obtained from equation (3.31): 

Using the explicit expressions (4.17) in equation (5.1), with = 1, we obtain 
to the following results: 

= :^b\r,)T, (5.5) 

We note that k± and k\\ are simply related by Z'^k,^\ = a.s it can be read off 
from equation(5.1), using equation (4. 15) and (4.1) and the fact that the small fre- 
quency behavior of the wave function is the same for both transverse and longitudinal 
directions. 

In the conformal limit b{r) = i/r, f{r) = 1 — (vrTr)^ and Tg = T/^/^, = 
1/ (nTy/^), we recover the results of [14, 15] for the holographic TV = 4 SYM: 

«:±Ar=4 = n^/X^a'^^T^ (5.7) 

where (i/is)^ = XN=i is identified with the A/" = 4 't Hooft coupling, in the AdS^ 
background. 

Finite Mass Quarks For massive quarks, the appropriate boundary condition 
should be Rirq) = 1, where rq is the UV cutoff determined by the value of the 
quark mass Mg, using equation (4.3). In this case, equation (5.3) gets modified and 
reads 



47rT. 



^R={^^^^] [l + Oiu)]. (5.9) 

Nevertheless, as in the conformal case [18], the results for k,± and remain 
unchanged, since they are independent of tq. Indeed rq enters in the wave function, 
as equation (5.9) shows, and cancels out in the w — )■ limit as we take the product 
in (4.15). 

5.2 The jet-quenching parameter 

As discussed in Section 2, the jet-quenching parameters can be defined in terms of 
the diffusion constants as: 

= 2 — , gil = — . (5.10) 

V V 
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The first parameter defines the transverse momentum broadening of heavy quark 
probes. 

There is also a different definition of the jet-quenching parameter, which is related 
to the perturbative relation between this quantity and an appropriate limit of a 
Wilson loop joining two light-like lines in Yang-Mills theory (see e.g. [47]). This was 
the basis of a different holographic calculation of g, that was carried out in [8] in the 
conformal case, and in [11] for the general backgrounds (3.1), which gives: 



where the subscript WL is introduced to distinguish this definition of the jet-quenching 
parameter from the original definition (5.10). As in the conformal case, this result 
differs from the result obtained via the Langevin equation. The reasons for this were 
analyzed in [53]. 

5.3 The diffusion constants via the membrane paradigm 

The method of the membrane paradigm that we introduced in section 4.2 allows us to 
read off the diffusion constants directly from the action (3.33) with no need to derive 
and solve for the fluctuation equations as in the previous section. The diffusion con- 
stants are defined in terms of the symmetric Green's functions by (5.1). Employing 
the relation (2.14) between the retarded and the symmetric Green's functions and 
the basic formula of the membrane paradigm, equation (4.22), we arrive at. 



where a = {_L, ||}. In the second equation above we used the fact that the metric 
dependence in (4.22) drops out in 2d. 

From the expressions (4.23), and using the near-horizon limit of the function Z 
from (5.4), we find, unsurprisingly, the same result as equation (5.5,5.6). 

We note that these results establish one of the very few examples of trivial flow as 
defined in [45], in the sense that the effective couplings Q, determined on the horizon 
membrane, stay unchanged through the flow from the horizon to the boundary. 
Therefore the dual field theory quantities i.e. the diffusion constants, which should 
be evaluated on the boundary, can also be computed directly on the horizon due 
to trivial flow. The other basic example of trivial flow is the shear viscosity rj/s. 
The reason for trivial flow is that there is no mass term in the fluctuation equations 
(3.36), because the geometry is fiat on the domain- wall directions. For the same 
reason one expects trivial flow for any transport coefficient that stem from string 
fluctuations on a generic domain wall background, as long as the fluctuations do not 
involve the radial direction. 




(5.11) 



(5.12) 
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5.4 A universal inequality: K| > k± 

From the expressions (5.5,5.6), one derives the ratio. 



„ , \,7 . , (5.13) 



«± \f'{rs)J f'{r,)h{rs 

where in the last equation, we used the definition of the world-sheet temperature T<j 
in (3.23). 

We note that the second term on the RHS of (5.13) is always positive definite in 
the deconfined phase T > T^. This can be seen as follows: First of all, f'{r)h{r) is 
a negative definite quantity at any r. This follows from the general relation, see e.g. 
[36]. 

/'W = -^M-)-'^ (5-14) 

p c 

where s is the entropy density and fe^; is the Einstein frame scale factor. The left 
hand side is manifestly negative definite. 

Secondly, the quantity h'{rs) is also negative-definite in the deconfined phase. 
This follows from the fact that, in the type of geometries that confines color, the 
string frame scale factor h{r) at zero temperature always possesses a minimum at 
some point r = r^,. Hence h'{r) < for r < r^, and h'{r) > for r > r^,. Moreover, 
we can argue that, for T > T^, the location of the bulk horizon should be closer 
to the boundary than r*, i.e. r/^ < r*, otherwise the Wilson loop would have linear 
behavior, as a result of saturation of the corresponding string at r^,. Since, the 
world-sheet horizon is always smaller than the bulk-horizon, it follows that, in the 
de-confined phase: 

rs<rh<r^, for T > T^. (5.15) 

Therefore b'{rs) should be negative-definite, and the entire second term on the RHS 
of (5.13) should be positive-definite. Therefore, we arrive at the universal result that 

K|l > K± for T>Tc (5.16) 

Equality is attained for f — 0. We check by a numerical computation in section 6 
that this inequality is obeyed in the particular background used in that section. 

5.5 A generalized Einstein relation 

We may derive a generalized Einstein relation by relating the diffusion constant 
(5.5,5.6) with the friction coefficient r]D in (2.17). 

On the one hand, we have found holographically the relation between the dif- 
fusion coefficients and the friction coefficients rj^^ for the Langevin equations in 
position space (2.24-2.25). From equation (4.18) and the definitions (2.13) we arrive 
at: 

k'^ = 2Ts7]'^ (5.17) 
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On the other hand, we can relate t]^^ to the momentum diffusion coefficients r/^, by 
equations (2.26). Therefore we find: 



= 2T,-fMgri^, (5.18) 

(5.19) 



p='-fMqV 



These relations lead to an important consistency condition: notice that, by equation 
(2.23), the coefficient 77!, must coincide with the zeroth order drag coefficient (3.18), 
calculated via the classical trailing string solution. 

Therefore, consistency requires that, inserting the expression (5.20) for 77!, in equation 
(5.19), the resulting expression for agree with equation (5.6). This is indeed the 
case: using the explicit expression (5.20), and the definition /(r^) = v^, we find: 



VD+P- 



dp 



^v \"'; , (5.21 



Inserting this expression in the right hand side of equation (5.19), and using the 
identity (5.13), the resulting expression exactly agrees with obtained from the 
Langevin correlator, (5.6). 

Finally, from equation (5.18) and the explicit expression (5.5), we find that rjj^ 
also equals the drag force coefficient (5.20). This implies that, at the end of the day, 
the friction term in the momentum diffusion equation is isotropic. 

The last equality defines the momentum diffusion time r = I/z/d. 
We arrive at the generalized Einstein relation: 

TK± = 2Mg-fTs. (5.23) 

where Mg is the quark mass and Tg is the emergent "world-sheet temperature" . (5.23) 
can be viewed as a generalization of the usual non-relativistic Einstein relation which 
has the form: 

TK = 2MgT. (5.24) 

The modified Einstein relations (5.18) and (5.19) had already been found in [20] in 
the particular case of the holographic dual of the M = 2* theory, and the results of 
that analysis fit consistently in the general framework we are considering here. 
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We note that in the conformal hmit, the world-sheet temperature is related to 
the bulk temperature as 

T 

Tg = , conformal, (5.25) 

and (5.23) becomes, 

Tn± = 2MgT^. (5.26) 

This is quite different from (5.24) and reduces to it only in the non-relativistic limit 
7 1. 

The generalized relation in (5.23) is defined in terms of a set of physical boundary 
quantities, and the geometric quantity Tg. In a sense, is the temperature that is 
read by the quark as it moves through the medium. provides the answer to 
the following interesting question: what is the temperature read by a thermometer 
moving with speed v inside a strongly coupled QGP of temperature T. 

Note that in the conformal case the answer is universal, and transforms simply 
with boosts. In the non-conformal case, associated with the Einstein-dilaton system 
the answer is dynamics dependent. 

5.6 Special limits of the diffusion constants 

In this section we study the diffusion constants (5.5) and (5.6) in the extreme rela- 
tivistic and non-relativistic limits and express these quantities in terms of thermo- 
dynamic functions. 

5.6.1 Non-relativistic limit 

As f{rs) = v"^, in the non-relativistic limit v 0, the world-sheet horizon approaches 
the bulk horizon: Ts — )■ r/j. Using the near- horizon expressions for the metric func- 
tions, one also finds from (3.23) that Ts — )■ T in this limit. Finally, we use the 
expression that relates the entropy density with the scale factor at the horizon, [36] , 

- = ^ = 4vrMX&lK), (5.27) 

to obtain 



2 /457r\ 5 £^ / s 



where Xh is the horizon value of A. On the other hand, the A/" = 4 result becomes 
i^±jv=4 vta/X^^T^. Hence the ratio becomes, 

^-A-r) ™ K, v^O. (5.29) 



A similar analysis for the parallel component yields 



2 /457r\ ^ £V s \ ^ _ 4 
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We note that this is exactly the same as (5.28). This is what one expects from the 
physical perspective. In the non-relativistic limit, the main source of momentum 
broadening is due to thermal fluctuations in the plasma, that itself is isotropic. 
Similarly the ratio of QCD and J\f = 4 results also become the same as in (5.29). 

5.6.2 Ultra-relativistic limit 

We consider the opposite limit f — )■ 1. Here the expression /(r^) = v"^ tells us that 
Ts —> hence the world-sheet horizon approaches the boundary. Using the near- 
boundary expression for /(r) in equation (3.6), and the near-horizon expressions for 
the metric functions, we find that, 

47rT, rsiiC/ey/^ (l + 0(log-i(r,))) . (5.31) 

where the constant C is given in equation (3.7). Upon substitution in (5.5), we finally 
obtain. 

We observe that the result diverges in the extreme relativistic limit v = 1. However, 
one obtains a finite expression by considering the ratio with the Af = A result: 

, (45vr^)f ^Y s \ H-'ilog{l-v')r oo^ 
Similarly, for the parallel component one finds. 



Again this is divergent as f — )■ 1, but the ratio with J\f=4 result again remains finite 
in this limit: 

(45vr^)f£Y . \^ (-|log(l-^^))-^ 

We observe that the parallel and perpendicular components of the diffusion con- 
stants asymptote essentially to the conformal result both in the f — > and in the 
f — )■ 1 limits, modulo logs in the second case and the appropriate adjustment of 
relevant parameters like temperature, entropy etc. 

We should warn the reader that, as f — )■ 1, we expect the break down of the treat- 
ment based on the local Langevin equation: in fact, for very large v the world-sheet 
temperature drops to zero, and the auto-correlation time of the Langevin Green's 
functions, Tc ~ T~^, diverges. On the other hand, the relaxation time td = tid—1 
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stays approximately constant. Therefore, eventually the relation td ^ r^, necessary 
for the local treatment, will break down at large enough v. In the next subsection 
we give a detailed discussion of this validity condition. 

Another important caveat, when considering the extreme relativistic limit of our 
results, lies in the fact that they are sensitive to the UV region of the geometry, 
and as discussed at length in previous work (see e.g. [37]), the details of the gravity 
theory we are using are not fully reliable in this limit. 

5.7 Time scales and validity of the local approximation 

The results of this section so far were obtained based on two separate approximations 
concerning the time scales involved. On the one hand, we assume we are in a short- 
time approximation, compared to the typical relaxation time. This means that, 
the quark velocity v can be assumed to be constant only within time scales that 
are much shorter than the relaxation time T£, = l/rjD. On the other hand, the 
analysis based on the local Langevin equation relies on a long-time approximation, 
this time compared to the typical time scales entering the Langevin correlators, and 
determined by the inverse temperature that the quark "feels" as it travels through 
the plasma. According to our previous discussions this is given hj Tc = 1/Ts where 
Tg is the world-sheet temperature (3.23). Therefore our analysis in terms of diffusion 
constants will be valid only for time scales t such that t^. t td. Existence of 
time intervals satisfying this condition requires that: 

11 . ^ 

— » 7f7- 5.36 

Since both rjD and T, depend non-trivially on the quark momentum, this condition 
translates in an upper bound on the quark momentum (or velocity), above which the 
local treatment breaks down^^. 

The relaxation time l/rjD is given by (3.18). We can write equation (5.36) more 
explicitly as 

We can read this condition as a lower bound on the quark mass. It is most 
restrictive in the UV (large f), where it reads, 

M,»^m^^m'. (5.38) 



^^Note that this condition is different than the condition given in [12] for the classical non- 
relativistic Langevin dynamics that is I/t? 3> 1/T. As the thermal behavior of the Green's functions 
are set by Tg rather than T, it should be this effective temperature that enters in the validity 
condition (5.36). 
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where C is defined in (3.7) and in the UV region A(rjj) is approximately given by, 

^'-) - - - Mog[A.V(C^)] - f^'^^' 

We observe that \{rs) vanishes in the extreme relativistic hmit because of the de- 
pendence on 7^^. However, this logarithmic dependence is milder when compared to 
the explicit dependence on 7 in (5.38). 

Alternatively, equation (5.36) can be read as a condition on the quark velocity, or 
momentum. For fixed quark mass, and for w — 1, approaches the AdS boundary, 
and the right hand side of equation (5.37) scales approximately as 7^/^. Explicitly 
we find 



v/^«2M,(^^J y^V^\j,) ■ (5.40) 

Consequently the condition (5.36) puts an upper bound to the quark velocity. 

To obtain an estimate of the upper bound on momentum, in the right hand 
side of equation (5.37) we can approximate the scale factor as 6(r) ~ A^/'^(r)(£/r)^, 
and replace the quantities and by by the corresponding conformal expressions, 
equations (3.24). We arrive at the bound (for an ultra-relativistic quark): 

For V close to unity, the dependence on v in \{r^ is very mild and the right hand 
side can be considered as a constant, depending only by the quark mass, temperature 
and value of the (holographic) coupling. We will give a numerical estimate of these 
quantities in the next section, in Improved Holographic QCD. 



6. Improved Holographic QCD and comparison with data 

In the previous sections we obtained general results for the correlators of world-sheet 
fluctuations, and for the Langevin diffusion constants, valid in any 5D Einstein- 
Dilaton theory admitting asymptotically AdS black hole solutions. Here, we will 
study in detail these results in Improved Holographic QCD, with the potential pro- 
posed in [38], whose thermodynamic properties are in good agreement with lattice 
YM thermodynamics [39] as well as the T = spectra of glueballs obtained on the 
lattice. 

We take the potential to be, as in [38, 11]: 

V^(A) = ^ {1 + K.A + V,\'l^ [log (1 + V,\'l^ + V^sA^)]'^'} . (6.1) 

^^In fact, the numerical studies in the next section shows that the other factor in the log dominates 
over the 7 dependence and A can be treated as a constant except in the extreme relativistic limit 

V = 1. 
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The coefficients Vi entering the potential, are fixed as follows (for a detailed 
discussion see [38, 11]): 

Vo = l^o. Vi = 14, ^2 = /3o^(^^^^fi^)', V^3 = 170. (6.2) 

where /3o and /3i are the first two pure Yang-Mills the beta-function coefficients, 

22 SI 

The coefficients Vq and V2 are fixed to match the perturbative YM /3-function, 
whereas Vi and V3 are fixed phenomenologically by comparing the equation of state 
of the model with that of YM on the lattice. 

The coefficient i is the scale of the asymptotic AdS^ space-time at r = and it 
sets the energy scale in the field theory. All observables defined holographically using 
the metric in the Einstein frame are measured in units of i. For a given class of black 
hole solutions with fixed UV asymptotics, the value of i can be set by matching the 
mass of the lowest glueball excitation [38]. 

It may seem that there is an extra scale associated to these models, with respect 
to 4D Yang-Mills, where there is a single scale, i.e. the quantity Aqcd setting 
the scale of conformal symmetry breaking in the UV. In our model, the analog 
of the QCD scale emerges as an integration constant that labels different solutions 
(distinguished by different UV boundary conditions) of the same theory, with i fixed. 
Exphcitly, it controls the UV asymptotics of the field A(r), given in equation (3.6), 
as A(r) = — (feologrA)"^ + 0(log~^rA). Therefore, it may appear we have two 
independent scales, i and A. However, as shown in [33], physical observables depend 
on A only via an overall scaling. Therefore, we can choose an arbitrary value for 
the dimensionless parameter £A, and subsequently fix the value of i to match some 
reference energy scale, as explained in the previous paragraph. 

The quantities that are computed by probe strings depend on another scale, 
independent of i, namely the fundamental string scale is- In string-derived models, 
the ratio i/ig is known. In phenomenological models on the other hand it must be 
adjusted to fit observation. For example, the ratio i/ig can be fixed by comparing 
the confining string tension of the holographic model (controlled by the Nambu-Goto 
action (3.8) , hence by ig) to the lattice value o"c = {AAOMeV)"^. One finds: 

6.5. (6.4) 
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This fixes the ubiquitous overall coefficient entering in the diffusion constants, in 
equations (5.5-5.6). 

With the potential given by (6.1), we numerically solve Einstein's equation for 
the metric and dilaton functions 6(r), /(r), A(r), to obtain black hole solutions of 
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(a) (b) 

Figure 1: (a) The ratio of the world-sheet temperature to the bulk black hole temperature, 
as a function of velocity, for different values of the bulk temperature. The dashed line 
indicates the ^dS'-Schwarzschild curve, = T/^/7. (b) The function ^Z(rs) as a function 
of velocity (with Z defined as in equation (3.31) and 7 = — f^), computed numerically 

varying the velocity, at different temperatures. The dashed line represents the conformal 
limit, in which 7Z = 1 exactly. 



different temperatures T, but obeying fixed UV boundary conditions (for a detailed 
discussion of the solution procedure, see Appendix A of [38]). Once the solutions are 
given, for each temperature we determine the position of the world-sheet horizon 
as a function of velocity v, by numerically solving the equation /(r.,) = v"^. 

From the metric coefficients evaluated at we can obtain the world-sheet tem- 
perature Tg, through equation (3.23). The ratio T^/T is plotted as a function of 
V, and for various bulk temperatures, in Figure 1 (a). We observe that Tg < T 
for all velocities, and as the bulk temperature increases this ratio approaches the 
AfiS'-Schwarzschild curve {Ts/T)Ads = (1 ~ v'^Y^'^. 

Another interesting quantity that provides an indication of how much the back- 
grounds deviate from the conformal case, is the function Z{r), defined in equation 
(3.31). For ArfS'-Schwarzschild, this function is exactly constant, Z{r) = 1/7. In 
Figure 1 (b) we portray the behavior of '-fZ{rs) as a function of velocity, for dif- 
ferent bulk temperatures. We observe that as f — > 0, 1 this quantity asymptotes 
to unity, as can also be seen analytically from equation (3.31) by taking the limits 
Ts — >■ 0, r<i — )■ r/j. Again, as the bulk temperature increases, we move closer to the 
AdS'-Schwarzschild behavior, represented by the dashed line in the graph. 

At this point, we note that Figure 1 also provides a confirmation of the universal 
inequality derived in section 5.4. The function Z{r) of equation (3.31) at the world- 
sheet horizon is given by Z{rs) = f'{rs)/4:'KTs. Then, from equation (5.13) we see 
that h\\/k± = Z{rs)~'^. On the other hand, the numerical computation shown in 
figure 1 implies that, 

Z{r,) <- Z{r,Y^ > 7' > 1. 

7 
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Therefore the numerics confirm that the inequahty Ky/K^ > 1 is satisfied. 

Knowing the numerical black- hole solutions and the values of r<., we can im- 
mediately compute the diffusion constants, (5.5-5.6). The results are discussed in 
subsections 6.2, 6.3. On the other hand, in order to compute the full Langevin cor- 
relators, we additionally need to solve the world-sheet fluctuation equations. This is 
analyzed in the next subsection. 

6.1 Correlators and spectral functions 

The retarded correlator of the trailing string fluctuations is given by equation (4.2) as 
a function of the frequency cu, where the wave-functions \E'^'"(r, oo) are the eigenmodes 
of the world-sheet fluctuations. From the full retarded propagator, one can further 
obtain the symmetric one through equation (4.18), and the spectral density through 
equation (4.19). 

To compute the Green's function (4.2) numerically, one must solve the linear 
fluctuation equations (3.37-3.38), with infalling conditions (5.2) at the world-sheet 
horizon r^, and unit normalization at the UV boundary. 

The numerical computation makes use of the shooting technique from the world- 
sheet horizon (specifying the in-falling initial conditions for the wave function and 
its derivative). Once the full solution is obtained, we normalize it dividing by the 
value of the solution at the boundary, in order for the wave function to obey the 
required boundary conditions. The results of the numerical analysis are shown in 
figures 2 through 4, and we will discuss them below in more detail. 

As discussed in Section 4, the real part of the correlator is UV-divergent and 
therefore it is very sensitive to the cut-off used in the numerical integration. Even if 
one subtracts the divergent term ~ 'yuj'^/e, one should be very careful to extract the 
limit e — )■ from the numerics, and eliminate all terms which grow as higher powers 
of u, but whose coefficient would vanish at e = 0. For example, after subtracting 
the divergence, the numerical calculation will be dominated by the subleading term 
in eq. (4.13), which at finite e grows as u^, but which is absent when the cut-off is 
removed. 

The imaginary part of the propagator on the other hand does not present these 
problems, and can be obtained in a clean way from the numerical computation. 
Since the imaginary part of the retarded correlator is related in a simple way to the 
symmetric correlator through equation (4.18), we chose to only show plots of the 
latter, which are discussed below. 

In Figure 2 and Figure 3 the symmetric correlator corresponding to a quark 
with infinite mass is shown, and it is compared to the WKB result in (4.40), for the 
transverse and longitudinal modes. From these plots we observe that the WKB result 
is a very good approximation to the spectral densities even at low frequency, which is 
a priori unexpected. In particular, from comparison with the WKB result, we learn 
that the symmetric correlators scale with a cubic power-law at large frequency. The 
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(a) (b) 

Figure 2: (a) The symmetric correlator of the _L modes by the numerical evaluation (solid 
line) and by the large- frequency WKB computation of Section 4.4 (dashed line), both in 
the Mq — 7- oo limit, (b) Difference of the numerical and WKB results. We show in each 
plot the curves corresponding to the velocities v = 0.1,0.9,0.99 and different plots for the 
temperatures T = Tc, 2Tc, 3Tc. 

difference w.r.t. the WKB result is small compared to the value of the correlator 
(the apparent discontinuity in some of the curves is an artifact due to the logarithmic 
scale — it is a small bump if plotted on a linear scale). 

Figure 4 shows the result for the finite mass correlators and their comparison to 



-45 - 




rGeV' 
^ fm 



10 15 20 25 30 



r = 3. T, 




10 15 20 25 30 



T=l.T, 



rGeV^ 
^ fm 




10 15 20 25 30 



Gsym"~GwKB 




10 15 20 25 30 



T = 3.T, 



-Gwkb" 



GeV^ 
L fm 




10 15 20 25 30 



(b) 



Figure 3: (a) The numerical result for the symmetric correlator of the || modes (solid 
line) together with the large frequency result from the WKB computation of Section 4.4 
(dashed line), for Mq = oo. (b) Difference between the exact and WKB results. As for the 
_L modes, we show in each plot, T = Tc, 2Tc, 3Tc, the curves corresponding to the velocities 
V = 0.1,0.9,0.99. 



the WKB approximation for large frequencies, using the results of section 4.4 and 
evaluating the factor (??). The correlators indeed display the linear behavior for 
large u that was derived analytically in section 4.4. 

The numeric computation in this case is performed by normalizing the wave 
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Figure 4: The retarded correlator ((a) real and (b) imaginary part) for finite but large 
quark mass, calculated numerically, (c) The symmetric correlator for finite quark mass, 
from numerical evaluation (solid lines) and WKB result of Section 4.4 (dashed lines), (d) 
The difference between the numerical and WKB result. In each plot we represent the 
curves corresponding to the velocities v = 0.1,0.9,0.99. We have taken T = 3Tc and 

Mq = Mcharm- 



function at the cutoff tq. The value of tq is determined by the quark mass (using 
equation (4.3)), rather than being given by the regulated boundary e. For the charm 
and bottom quarks, we take Mcharm = 1.5 GeV and Mbottom = 4.5 GeV.^^ 

The correlators are also evaluated at rg, following the formula (4.2). We use the 
results obtained in [11] for the cutoff, yielding tq ~ 1.4 for the bottom quark and 
rq ~ 7.5 for the charm quark. We chose to show in Figure 4, as an example, the 
results for the transverse mode of a charm quark at T = 3Tc, for different velocities. 



^^These values are subject to renormalization in the plasma due to interactions with the media. 
Therefore, they become temperature dependent. However, as we show in [11], this temperature 
dependence is very mild in our holographic model, within the relevant temperature range. 
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Figure 5: The ratio of the diffusion coefficients k± and Ky to the corresponding value 
in the holographic conformal = 4 theory (with A_/\/=4 = 5.5) are plotted as a function 
of the velocity v (in logarithmic horizontal scale) from equations (5.5)-(5.6). The results 
are evaluated at different temperatures T = Tc, 1.5Tc,3Tc in the deconfined phase of the 
non-conformal model. 



6.2 The jet-quenching parameters 

The diffusion constants k± and /ty are computed directly, as a function of temperature 
and velocity, by evaluating equations (5.5-5.6) at the world-sheet horizon specified 
by equation (3.15). 

To give an idea of the effect of the running of the scalar field, and of the breaking 
of conformal invariance in our model, in Figure 5 we show the ratios of k± and k\\ to 
the corresponding quantities obtained in the AdS black hole background, equations 
(5.7-5.8) representing strongly coupled A/'=4 SYM. The ratios n/nconf are shown as a 
function of velocity, at different temperatures. The conformal results are obtained by 
fixing is by its AdS /CFT relation to the fixed coupling of AT = 4 SYM, (i/is)^ = A. 
We take as A = 5.5 as in [5]. 

From Figure 5.1 we observe that, apart from an overall normalization, the non- 
conformality in this particular model significantly affects the diffusion constants only 
for temperatures close to T^, and for velocities that are not too large. Indeed, if we 
choose A ~ 0.5 instead of A = 5.5 in the conformal case (just to make the overall 
magnitudes similar in the comparison), then our result agrees with the conformal 
result within the 10% level, in the range v > 0.6 and for T > 1.5Tc. 

In the rest of this section we will focus on the jet-quenching parameters q, which 
are related to the diffusion constants by 

.1 II f^^^ , \ 

q^ = 2 — , gll = - . 6.5 

V V 

The Langevin dynamics defines the two independent parameters, q± and gy. The 
first controls the transverse momentum broadening of a heavy quark probe moving 
through the plasma, and it is the one usually quoted in relation to experimental 
results. In this Subsection we present the result for the dimensionless quantities q/T^ 
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Figure 6: The jet-quenching parameters g_L and obtained from the diffusion constants 
(5.5)-(5.6), normalized to the critical temperature Tc, are plotted as a function of the 
velocity v (in a logarithmic horizontal scale). The results are evaluated at different tem- 
peratures. 

since they do not depend on how we fix the overall energy scale in the holographic 
QCD model. We will translate the result to physical units in the next Subsection. 

In figure 6 we plot the two jet-quenching parameters (in units of the critical 
temperature T^) as a function of the velocity, for different values of the temperature 
T. The behavior for small v is dominated by the 1/v factor in the definitions (6.5). 

We note that the difference between the longitudinal and transverse modes is due 
to the function Z defined in equation (3.31), which reduces to f'{rs)/4:7iTs when eval- 
uated at the world-sheet horizon. More precisely, the relation between the diffusion 
constants is k,± = K,\\Z'^{rs). In sections 5.6.1 and 5.6.2 we showed that 'yZir^ — 1, 
both when — and v ^ \. This is also apparent in the numerical result shown in 
Figure 1 (b). 

It is instructive to translate the quark velocity v on the horizontal axis of Fig. 6 
into momentum, p = M^f 7, where Mg is the quark mass. Taking Mcharm = 1-5 GeV, 
Mbottom = 4:.5GeV, the resulting plots are shown in Fig. 7^"^. From these plots, we 
observe that q± is almost constant over a wide momentum range, for a relativistic 
heavy quark. This is not so for gy, the difference being due to the extra factor of 
Z'^i^s) ~ in the latter. 

From Figure 7 we observe that, for a fixed momentum, q increases with tem- 
perature, as can also be inferred from the analytic expressions (5.5)-(5.6). This 
behavior is shown more clearly in Fig. 8, that displays g as a function of tempera- 
ture (in units of the critical temperature Tc) for different quark momenta. This is 
the q± oc behavior predicted by both the relativistic approximation (5.32) and 
the non-relativistic one (5.28), once we use the fact that s ocT^ approximately. 

It is important to keep in mind that, for a finite quark mass, the trailing string 
description breaks down when the world-sheet horizon coordinate becomes smaller 
than rq, the point where the trailing string is attached (this is infinite for an infinite 

^'^As shown in [11], thermal corrections to its mass are negligible in our set-up. 
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Figure 7: The quantities q_\_/T^ and (^y/T^ plotted as a function of the quark momentum 
p. The plots for the charm and the bottom quark differ by a scaling of the horizontal 
direction. 



quark mass). This translates into an upper bound on the velocity, or momentum, 
beyond which the setup breaks down, and the string dynamics on the flavor branes 
becomes important. To estimate this bound, we can use the AdS^ relations for 
the quark mass and world-sheet horizon: for a heavy quark rg is approximately 
tq ~ (£^/27r£g)A^/^(rQ) Mg^, and for large momentum p, ~ 7~-^/^(p)(7rT)^-'^, with 

garania{p) = ^1 + p'^/M^. This results in the approximate bound: 

Notice that the bound is stronger for higher temperatures and lower quark masses. 

We estimated numerically the bound on p in the model we are using. The results 
are displayed in Figure 9. From this figure we see that the the bound is easily satisfied 
for both the Charm and Bottom quarks in the RHIC and LHC regimes, if we use 
Tc ~ 200 MeV 

Finally, we need to check what is the allowed range of p such that the local 
approximation to the Langevin equation is reliable, as discussed in Section 5.7. For 
this to be the case, we need the quantity Ts/rjD to be large. For an ultra-relativistic 
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Figure 8: The jet-quenching parameters q±_ and q\\ plotted as a function of T/Tc, for 
different quark momenta. 
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Figure 9: This figure displays the upper bound on momentum Pmax (in logarithmic scale) 
beyond which the trailing string picture ceases to be valid, as a function of the quark mass, 
and for different temperatures. 



quark, this condition translates to equation (6.6). In our numerical solution we obtain 
As ~ 3 X 10~^. Therefore we expect that, for moderate temperatures the bound is 



pretty mild^^. Taking the above result as a reference value for A^, and ^s/^ — 0.15 
from equation (6.4), we can rewrite the bound (6.6) more explicitly as follows: 



For example, for the charm quark (M^ ~ 1.5 GeV), and close to the critical temper- 
ature Tc, this translates into p <^ 2 Gel^(1.5 GeV/Tc)^. 

However, the situation changes dramatically as temperature increases: from eq. 
(6.7) we observe that the upper limit at temperature T decreases as (Tc/T)^. A 
graphical representation of the validity condition is shown in Figure 10, computed nu- 
merically from equation (5.41) for both the charm and bottom quarks. The p- region 
in which the diffusion process can be approximated by a local Langevin equation, 
with constant friction and diffusion coefficients, lies in the left side of the vertical lines 
(each corresponding to a different temperature). From these plots we observe that 
the bound is satisfied for momenta up to ~ 70 GeV (charm) and more than 200 GeV 
(bottom) at T = T^, but for larger temperatures the bounds are much stronger. For 
example the bounds at T = 3Tc are Pcharm < 10 GeV and Pbottom < 100 GeV. 

What these values of T correspond to in terms of actual physical temperature of 
the QGP, is a subtle question, as we will discuss more in detail in the next subsection. 
However an order of magnitude estimate can be obtained by setting ~ 180 GeV 
in these plots. This means that, for RHIC temperatures and momenta, the local 
approximation remains valid. On the other hand, if the holographic setup is to be 
applied to ALICE results, it is likely that one should use the full non-local form 
of the generalized Langevin equation, and the simple parametrization of transverse 
momentum broadening in terms of q breaks down. 

6.3 Comparison with heavy-ion collision observables 

Fit of RHIC data for nuclear modification factors with hydrodynamic simulations 
prefer a strong jet-quenching parameter for light quarks about q± = 5 -15 GeV'^/fm 
(for a review of recent results and a more references, see e.g. [47]). 

In order to compare our results to QGP observables we need to evaluate the 
results of the previous section at typical temperature for QGP Tqgp ~ 250 MeV. 

However, as discussed in detail in [11] it is not easy to make a direct comparison, 
because our calculations are made with a pure glue background (neglecting therefore 
the quark contributions). It was argued recently that quarks contribute importantly 
in energy loss, beyond their enhancement of the number of degrees of freedom, this 

^^There is a certain degree of arbitrariness in the choice of normahzation of A. However, changing 
the normahzation of A would result in a value of £/£s different from the one we are using here, 
equation (6.4). The important thing is that, ones we insist in fixing the confining string scale at a 
certain physical value, the quantity £/£s scales as X~^^^ under an overall scaling of A [38]. Therefore 
the bound (6.6) is independent on the overall normalization of A(r). 




(6.7) 
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Figure 10: The quantity Ts/rju is plotted against quark momentum, for different bulk 
temperatures. Figures (a) and (b) refer to the charm and bottom quark, respectively. For 
each temperature, the validity of the local Langevin equation constrains p to the left of the 
corresponding vertical line, which marks the transition of Ts/rjo across unity. 



was shown to be the case in the thermal D3 — D7 system [54], and the same feature 
was already noted in [56] in a non-critical model and a in a model based on wrapped 
D5 branes. 

To proceed further we will translate the physical QGP temperature to our T. 
To do this requires picking up a comparison scheme. In a direct scheme one simply 
takes T = Tqgp- 

On the other hand, one can argue that the relation between the QGP temper- 
ature and that of the holographic model should be such that the energy densities 
are the same. Energy density scales as the number of degrees of freedom, and the 
holographic setup we study is supposed to describe pure Yang Mills theory, rather 
than QCD with three light flavors. Therefore, matching energy densities leads to 
a holographic temperature T higher than the QGP temperature, due to the differ- 
ent number of degrees of freedom in the two theories. This reasoning leads to the 
identification of an alternative scheme, referred to as the energy scheme^^ , where the 
effective temperature T is related to the real QGP temperature Tqgp by the implicit 
relation [11]: 

^hol{Tenergy) — 11-2 Tg^p (6.8) 

where ehoi{T) is the energy density of the holographic model. 

Computing the ehoi{T) numerically one obtains the approximate linear relation 
(as shown in Fig 11 ) between the direct scheme and energy scheme temperatures: 



T 

-'■ energy 

MeV 



23.7+1.2 



dir 



MeV 



(6.9) 



^^One can also define an entropy scheme, where one matches entropy density rather than energy 
density. We checked that the numerical results obtained in the energy and the entropy schemes are 
essentially very close to each other. 
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Figure 11: Relation between the temperature in the direct and energy schemes. 
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Figure 12: The jet-quenching parameter q± in the direct (dashed hues) and energy (sohd 
hnes) schemes, as a function of temperature, for different quark velocities. 

The same Tqgp corresponds to a higher temperature of the holographic model 
in the energy scheme, than in the direct scheme. Therefore, using the energy scheme 
to match the QGP temperature results in higher values for q, than those obtained 
in the direct scheme. This behavior is apparent in Figure 12. 

We are now ready to translate in physical temperatures the results for q± pre- 
sented in the previous Subsection. This is done in Figures 13 and 14, which are 
analogous to Figures 7 and 8, except that the temperature and gl are displayed in 
physical units. In order to express q and T in GeV'^/fm and MeV, respectively, we 
have to introduce physical energy units. The overall energy scale was fixed, as briefly 
explained at the beginning of Section 6 (and in more detail in [38]), by matching one 
dimensionfuU quantity to its physical value. As in [38], for this purpose we used the 
lattice value of the lowest glueball mass. 

Figure 13 shows q± in both schemes, as a function of the probe quark momentum, 
for charm and bottom quarks and for various temperatures in the range relevant for 
RHIC and for the ALICE experiment at LHC We observe that, although the values in 
the energy scheme are higher, at the temperatures relevant for RHIC (T ^ 250 MeV), 
q± varies in the 1 — 5 GeV^/ fm range except at low momenta where it is substantially 
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Figure 13: The jet-quenching parameter q± in the direct (dashed lines) and alternative 
(solid lines) schemes, as a function of momentum, for different physical values of tempera- 
ture. 

higher. A full Langevin fit is necessary in order to ascertain if these numbers fit the 
data. 

For the highest temperature shown in Figure 13, T ~ 800 MeV (which is not in 
the range of RHIC, but may be within the reach of ALICE), the predicted value in 
the energy scheme reaches g ~ 25 GeV"^ / fm. 

We notice from Figure 13 that up to very large momenta q± is effectively in- 
dependent of p. Therefore one can safely neglect the non-linearity in the Langevin 
equation for large ranges of p. This also allows to pick a reference momentum (say, 
p PS 10 GeV) and study more closely the behavior of q± as a function of temperature. 
This is done in Figure 14, which is analogous to Figure 8 but with physical units for 
the temperature. 

In figures 13 and 14 one must keep in mind that the parametrization of transverse 
momentum broadening by q± fails for momenta that are out of the range validity of 
the long time approximation to the diffusion process. These bounds are given at the 
end of the previous section as a function of T/Tc. In the direct scheme, a plasma 
temperature of 250 MeV corresponds to T ~ Tc in the holographic model^^, but in 
the energy scheme the same plasma temperature corresponds to T ~ 325 MeV = 
1.3Tc according to the relation (6.8). Similarly, the temperature Tqgp = 800 MeV 
corresponds to T ~ 980 MeV ^ 4Tc in the energy scheme. Therefore, comparing 
with Figure 10, we deduce that at least for the charm quark, we cannot really trust 
the analysis in terms of q± for the highest temperature displayed in Figures 13-14, 
unless we look at momenta smaller than 4 GeV. On the other hand, the results 
shown for the bottom quark are consistent with the long-time approximation. 

^''We remind the reader that the critical temperature of the IHQCD model we are using is roughly 
Tc = 247 MeV [38]. 
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Figure 14: The jet-quenching parameter q± in the direct (dashed hnes) and energy (sohd 
hnes) schemes, as a function of temperature, for probe charm and bottom quarks with 
Pt ~ lOGeV. The bottom plots on show the same functions on a narrower temperature 
range than the top plots. 



Finally, we should remind the reader that our results in the UV are not totally 
trustable due to the fact that the background used may not describe properly the 
details of QCD. 
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Appendix 



A. Boundary terms and scheme dependence in the propaga- 



The form of the retarded propagator we have obtained in Section 4, equation (4.2), 
depends on the form of the action we took as a starting point, namely equation (3.8). 
As the correlators have UV divergences there is potential scheme dependence in their 
calculation that we now address. The associated counterterms do not modify the 
wave equations (3.36), but they change the value of the on-shell action and therefore 
the correlator. 

Since the boundary terms that we add to the action must be real, the scheme 
dependence can only manifest itself in the definition of IieGii{u!). Therefore, it does 
not affect the physical quantities described in Sections 4 and 5, such as the diffusion 
constants and spectral densities. In the case at hand, as we show below, the only 
scheme dependence in lie Gfi{uj) amounts to a renormalization in the (heavy) quark 
mass. 

In what follows we study the divergence structure of the action (3.8), expanded 
to quadratic order in the fluctuations defined in equation (3.26), around the classical 
trailing string solution. To regulate the action, we cut-off the r-integration at r = 
e > 0, and we study the divergences in the e — ?■ limit. 

Starting from equation (3.27) and (3.28), we obtain to quadratic order in the 
fluctuations. 



fore the integral is dominated by the region around r ~ e, and it is given 
approximately by: 



tor 



Sng — 5'o + 5*1 + 5*2 + ■ ■ ■ . 
Below, we write explicitly and discuss each term separately. 

• Zeroth order. 

The zeroth order term reads simply: 



(A.l) 





For small r, the integrand is approximately equal to (£/r)^ A^/'^(r)7 ^. There- 




(A.3) 



giving a 1 /e divergence. 



• First order. 

The first order term in the fluctuations is a boundary term of the form: 




(A.4) 
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Since on-shell 5X^^{r = 0,t) is finite, this term is not divergent. 



Second order. 

At quadratic order the action is given by equation (3.34), 



So 



1 f 1 

/ drdr -H""^ 



1 ^ 



with 



f)4 



i=2 







a/(/-^^)(6^/-C2) 

Vif-^'Wf-c^) 



(A.5) 



(A.6) 



It is convenient to write the fluctuation in Fourier space. The solution of the 
field equations, equations (3.36) reads, close to the boundary: 



(A.7) 



Inserting the above expression into the action, and using the asymptotic ex- 
pressions. 



- - l/7^ (Ai 



we observe that the only divergent term as e — ?■ originates from the terms 
involving two time derivatives of 6X: 



e 7 2 
_ A^/=^(e) f 1 
€72 



duco^j'\C^{u)\' + j'\C!{u;)\'' 



dt-i'^ [5X 



7 



4/^^11 



(A.9) 



We can reabsorb both divergences, (A. 3) and (A.9), with a single covariant 
boundary counterterm. 



S. 



count 



AM(e) j dt \/x^X^, 



(A.IO) 



which corresponds to a renormalization of the quark mass. Indeed, expanding equa- 
tion (A.IO) to second order in X = vt + 6X, we find: 



S, 



count 



AM(e) 
7 



-f^6X^)\-f^ (6X ^^ 



(A.ll) 



Comparing with equations (A. 3) and (A.9) it is clear that the following choice of the 
leading divergence of AM cancels both the leading and second order divergences: 



e 27r£2 • 



(A.12) 
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From the discussion above, we conclude that the only boundary term that can 
lead to a scheme dependence of ReG/j(a;) can come from the finite part of the 
counterterm (A. 10), with a finite coefficient 5m completely specified by fixing the 
renormalized quark mass. 

An independent way to see the same effect, is as follows. According to the first 
line in equation (A. 9), a finite counterterm of the form (A. 10) with coefficient 5m 
would shift ReGij(u;) by a term proportional to uo^: 

£2 £2 

AReG^H = — -5m7w^ AReG^^u) = —— 5m 'j^ u'^ (A.13) 
znij 2n£j 

In the Langevin equation, this amounts simply to a finite shift in the quark mass: 
the generalized Langevin equations (2.2) read, in Fourier space, 

-fMq5X^{uj) + Gji{uj)X^{uj) + ^^{uj) = 0, (A.14) 

a;2 7^M,(5Xll(cj) + g5j(cj)xII(cj) +^l'(a;) = 0. (A.15) 

after expanding to first order in fluctuations. 

The conclusion is that the shifts (A.13) are equivalent to a finite renormalization 
of the quark mass, Mq — >■ Mg + i'^/ {27ri1)5m. Therefore, once the renormalized quark 
mass is fixed e.g. at zero-temperature by fixing the counterterm, there are no further 
ambiguities in the two-point function. 



B. Analytic calculation of the diffusion constants 

Here, we provide a derivation of equations (5.5) and (5.6). These equations follow 
from the flux (4.15) which can be evaluated at any point, in particular at the horizon. 
Let us define, 

for notational convenience. 

Near the horizon, the solution of the fluctuation equations (3.37) and (3.38) are 
of the form: 

5X^ = Cx(r, - r)-^^, 5X\\ = ^.(r, - r)"^"^. (B.2) 

In (4.15) we also need the near-horizon expression for the r-r component of the world- 
sheet metric if''^ equation (3.35). Using the definitions C = vh{rs) and /(r^) = v'^ 
(3.20) and the definition of Tg in (3.23), we find: 

H"-' ^A'KT,h\r;){r-rs). (B.3) 

Substituting (B.2) and (B.3) in (4.15) yields the flux near the horizon (and every- 
where) : 

f \C^\H\r,)u, , ^ 

ImGi? = <^ 2 h\r,) , , (B.4) 
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From (3.31) we also find, 



f'(r ) 

^W^^, r^r, (B.5) 



Therefore tlie calculation is reduced to finding the coefficients C± and Cy . This can 
be done by matching the near-horizon solution (B.2) for small u to the exact analytic 
solution of the fluctuation equations again for small u. We give details for 6X±, the 
other component is entirely analogous. For small u (B.2) expands as, 

6X± ^C±- tooC± log(r, - r). (B.6) 

In the strict a; = limit this gives 

6X^ = Cx, 00 = 0. (B.7) 

On the other hand, (3.37) can be solved exactly in the strict a; = limit: 

5X^ = Ci + C2 / . (B.8) 

Jo Vif-v^Wf-C^) 

Requiring unit norm on the boundary fixes Ci = 1. The second term diverges at the 
horizon, therefore in the strict u = limit, one should have C2 = 0. Therefore one 
has, 

SX^ = 1, (B.9) 

in the strict u = limit. Since the solution (B.9) is valid everywhere, including the 
horizon, one can match it with (B.7) and obtain 

C± = 1. (B.IO) 

More generally, the solution at all r can be written as 

5X^ = C^{u) (r, - [l + Z^ix(c^)(r, - r) + D2±(c^)(r. - rf + O ((r, - rf)] . 

(B.ll) 

Expanding C± and D± around uj = (the regular expansion is guaranteed by the 
regularity of the a; = solution), we get 



SX^ = Cx(0) (r,-r) 



1 + (A±(0)(r, -r) + ((r, - r)^)) + 



(C;(0)/Cx(0) + D[^{0){r, -r) + ((r, - r)')) u + O [u^) 



,(B.12) 
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Now, equation (B.9) implies that -Di_l(0) = — and so on for all the Di±{0) — and 
C± = C±{0) = 1. Hence the solution for all values of the radial coordinate is (5.3) : 



SX< = (r, - r) 



(B.13) 



with Ci±{r) ^ C^(0)/C^(0) + D[^{0){rs - r), close to the horizon. 

In passing, we note that slightly away from the cu = limit one can allow for the 
second term in (B.8), and expanding the integrand near the horizon, one obtains. 

Matching this with (B.6) one can also determine C2 in the small u limit: C2 = 
—iujlP'{rs). This information is not required to calculate the diffusion coefficients. 

Use of (B.IO) in (B.4) and eventually in (5.1) yields the desired result (5.5) (after 
including the string tension in front of the world-sheet action). The discussion for 
the parallel component is similar. Solving the fluctuation equation (3.38) for 07 = 
and by matching (B.2) one finds C\\ = 1 and using this and (B.5) in (B.4) yields 
(5.6). 



C. Details of the WKB approximation 

We follow the steps outlined in section 4.4. It is convenient to define the dimensionless 
variables x = r/vg G (0, 1) and Us = ojTs. The analog Schrodinger equation is 

-r + v,{x)i, = 0, Vs{x) = -j^ + ^(log7^)" + ^(log7^)'^ (c.i) 

where 

and the functions R{x) and Z{x) are: 

R=Vif-v'Wf-C'), Z = bW{f-v')/ib'f-C^) (C.3) 

We divide the range < x < 1 in three regions, in each of which we use different 
approximations to solve the Schrodinger equation. The following discussion holds for 
both ± and || fluctuations, so we will not make any distinction from now on. 

1. Near Boundary: a; ^ 1 

In this region we have the following asymptotics: 

R{x) ~ 6^(x)/7, Z ~ 1/7, X < 1, (C.4) 
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and the Schrodinger potential is approximately 

K ^ -7V + (log&)" + (log^)'^ x^O. (C.5) 

In the near-boundary region the Einstein frame scale factor becomes that of 
AdS space-time, and we have: 

hix) ~ (C.6) 

One important property of this region, is that the quantity rA'/A is small. The 
reason is that the field equation for A(r) is [32, 33]: 

A'(r)~^A2~- ^ r^~A«l (C.7) 
I r A 

where 6£;(r) ~ £/r is the Einstein frame scale factor close to the boundary. 

As a consequence, all terms proportional to rA'/A , or corresponding higher 
derivative terms in A, can be treated, to a first approximation, as subleading 
in an expansion in A. 

2. Near Horizon: x ~ 1 

In this region we have: 

R{x) ~ (47rT,r,) h\r,) (1 - x), Z ~ (^-8) 



47rT, 



V,{x)--[u^ + \\—^^, x^l (C.9) 



leading to: 

T/Jt"! ~ — I (~? -I- 

47 (1 - 



where u = u/inTs. 

3. WKB Region: xtp <^ x < 1 

This is the classically allowed region, where Vs{x) < 0. For large Ug, the first 
term in equation (C.l) dominates, except close to the turning point Xtp, where 
the contributions of the other terms get large, 

xtp<^x<l (C.IO) 

Since, for any large but finite Ug, Vs{x = 0) = +00, the turning point for large 
Us is close to the boundary and it is found by solving the equation Vs{x) = 
in this limit, i.e. using V in the form (C.5). Keeping in mind that derivatives 
of X{x) close to the boundary produce corrections of 0(A) ^ 1, we find the 
turning point for large Us'- 

Pi 

Xip = — (1 + 0(A)), u;, >1. (C.ll) 
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The crucial fact is that, for large Ug, Xtp -C 1, and regions 1 and 3 overlap. On 
the other hand, regions 2 and 3 overlap close to x ~ 1. Therefore, the solution 
in the WKB region can be used to connect the near-boundary and near-horizon 
asymptotics. 

To find the wave-function in the large ujs regime, we follow the steps outlined in 
Section 4.4. 

1. Consider first the WKB region. The two independent solutions to —ip + Vil) = 
in the region <^ are written, in the WKB approximation: 



ipi^ — cos / p, 



^2 ~ ^ sin j p, 



Explicitly, the general solution has the form: 



wkb 



R 



cos 



L^o R 



C2 sin 



R 
(C.12) 



xtp<^x<l (C.13) 



In the equation above, we made the arbitrary choice x = for the lower inte- 
gration limit, in order to avoid ambiguities in the definitions of the integration 
constants. 



2. Consider now the near-horizon region. There, we can solve Schrodinger's equa- 
tion with the potential (C.9). The solution with in-falling boundary condition 
at the horizon is 



iph ^ Ch{l - x) 



X] ^"^+5. 



UJ 



UJ 



X 



1, 



(C.14) 



Since both forms (C.13) and (C.14) are valid in the near-horizon region andioT 
large u, we can relate the coefficients by evaluating (C.13) near the horizon. In 
order to use the near-horizon expansion of the equation (C.8) in the integrands 
appearing in equation (C.13), we change the extremum xq to another point 
Xi belonging to the horizon region. This introduces a common phase shift, 
9 = J^^ uJsb'^/R, in the sine and cosine functions. Taking this into account, the 
near-horizon expansion of equation (C.13) reads: 



wkb 



{ATiT.rsY'^ (1 - xy/^{Ci cos[e - wlog(l - x)] + 
-hC2 sin[6' - a;log(l - x)]}, x 1. 



Comparing equations (C.15) and (C.14) gives the relations: 



-1C2 = Ci 



-id 



(47rT,r,)V2 



(C.15) 



(C.16) 
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Next, we consider the boundary region, a; <^ 1. Here the potential has the form 
(C.5). Since the potential diverges as a WKB treatment is impossible 

all the way to a; = 0, so we must resort to another method. The strategy we 
follow is that of an expansion in the derivatives of A(x), more precisely in the 
small quantity rX'/X ~ 0(A) ^ 1. This will allow us to write an approximate 
expression for the solution, valid for any Us- 

Using the approximation for the metric in (C.5), the second entering the near- 
boundary potential can be written as: 

(logfe)" ^ 1 + - o ^ h (log&) 




3 \ \ J X \ 3A 



(C.17) 



The A-dependent terms in the parentheses are 0(A) or O(A^). 

One may naively think that it suffices to solve Schrodinger's equation keeping 
only the leading terms in (C.17), and neglecting the 0(A) corrections. However, 
as we show below, these subleading terms affect the leading term in the solution. 

Let us ignore for the moment the terms containing derivatives of A in equation 
(C.17). them, the near-boundary Schrodinger equation with the potential (C.5) 
reads: 

-V^" + -V^ = (C.18) 



X 



whose general solution is: 

sm{'yusx) 



^° (x) = 



cos(7a;sx) 

sm['yUsX) + 



+ A2 



cos{'yuisx) 



'JUX 



(C.19) 



However, this cannot be the full story, for the following reason. For both small 
X and small 'yUsX, we can expand this solution as: 

ruvix) ^(—)--( ^u.x « 1 (C.20) 



^"yoosj X \ 3 

On the other hand, in the same regime 'jUgX — ?■ 0, we can ignore the constant 
term (7^;^)^ in the potential, and Scrhodinger equation becomes: 

< = (C.21) 

which has the exact solution: 



ij]jy{x) = Cs h{x) + O. h{x) ^ (C.22) 



where Cg and 0„ are integration constants corresponding to normalizable and 
non-normalizable solutions. Using the near-boundary form of the metric (C.6) 
this expression becomes: 

^\jy{x) ~ —\^''\x) + 0„ x^\-^'^{x), -fUsX < 1 (C.23) 
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and it does not agree with the small jUgX expansion (C.20) due to the extra 
factors of A^^^^. From this discussion, we conclude that to find the correct 
behavior of the solution near the boundary, one cannot completely ignore the 
terms containing A' and A" in equation (C.17). Therefore, the true solution in 
the boundary region, rather than (C.19), will read instead: 



source vev 

COfiipfUJsX) 



+A2 



sin(7a;sx) + 

sin(7Co'sx) 



COs(7CJsX) 



Fi(x,7c^.) (C.24) 



where Fi and F2 are some unknown functions, with asymptotics: 

F^{x-iUJs) ~ \^'\x), F^ix-fu,) ~ A-2/3(a;), ^uJsX < 1. (C.25) 

On the other hand, we know that the functions Fi and F2 must be replaced by 
constants in the limit when A(x) is not changing at all. This suggests that we 
can parametrize the functions Fi and F2 as: 

70;,) = X^^/^{x) [1 + ^i{x, 7W,)] (C.26) 

where the functions ipi^2{x,'yuJs) are small compared to unity for small xX'/X 
and the + sign in the exponent of A corresponds to z = 1, and the — sign to 
i = 2. 

One could in principle derive differential equations for (pi^2, and solve them per- 
turbatively. This is equivalent solving the fluctuation equation h" + (log R)' h' + 
ah = with some book-keeping parameter a and where the prime denotes 
derivative w.r.t. the variable x'yUs- The solution can be found perturbatively 
in a series expansion both in a and x, such that X{x) can be regarded as a 
small expansion parameter. A simple limit of this solution is to keep only the 
leading term in the latter expansion and sum up the perturbative series in a 
fully and then set a = 1 to recover the original fluctuation equation. This can 
be achieved in an iterative manner and the answer is indeed given by (C.24) 
and (C.25). This method justifies the appearance of the extra factors Fi and 
F2 in (C.24). 

However, proceeding with this method and combining it with the WKB ap- 
proximation in order to achieve the full WKB solution requires going beyond 
the leading term in the perturbative series in the expansion in A(x), which is 
very cumbersome. Therefore, in the following we will follow a different strategy, 
which will give us a simpler (albeit more crude) way to estimate the coefficients 
of the WKB wave functions. 
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4. In the limit ^ 1, the UV region x -C 1 overlaps with the WKB region 
X > Xtp , because x^p <C 1. Therefore, the UV solution (C.24) must match, for 
large '~fUJs, the small-x limit of the WKB solution (C.13), which using the UV 
expansion of R{x), equation (C.4), reads: 

Ci C9 
ipwkb - cos(7a;sa;) + — ^ •mvijbJsX), x < 1. (C.27) 
a/7 V7 

In order for (C.24) to agree with this expression, it is necessary that the func- 
tions F\[x) and -F2(x) become approximately constant for x ^ Xjp, where 
(C.27) can be trusted. Since Xt^ ~ 1/Ws7, we conclude that: 

Fi{x) -f Fi = const 'joJsX > 1. (C.28) 

z = 1, 2. In order to complete the matching we must estimate these constants 
Fi and F2. One way to argue is as follows: to satisfy equation (C.28), we need 
the functions ipi and defined in (C.26) to have the following property: 

Fi 

(fi{x) ^ -1 + ^±2/3(3;-) ' x>xtp (C.29) 

again, where the + is for i = 1, while the — for i = 2. But as we argued that 
the functions (pi and (p2 must stay small for slowly varying A, setting ifi^^O 
gives us an estimate: 

Fi ~ A2/3(a;o) , F2 ^ A-2/3(a;o) , x » x^p (C.30) 

where xq is a point in the vicinity of the turning point Xtp. Since A is slowly 
varying, we can take xq = Xtp , and the error we make will be of the order 
XtpX'/l ~ 0(A) < 1. 

Therefore, we match (C.27) with the large Ug limit of (C.24), in which the 
functions Fi and F2 are replaced by the constants A^^/^(xtp): 

0/0 r) In 

ipuv AiKp sm{jUsx) + A2\tp' cos{-fUsx) (C.31) 
where we have defined Xtp = \{xtp). For large u the turning point is given by: 

Xtp ^ , rtp ~ — . (C.32) 

Matching (C.31) and (C.27) we obtain: 

Ci = A,7A2v/7, C2 = \lpA^^, (C.33) 
which through equation (C.16) relates Ch to A\. 
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5. Finally, we determine Ai, and consequently all other coefficients, by imposing 
unit normalization of the function '^{x) = lZ~^/'^'il){x) aX x = ri,/rs^ i.e. the 
point where the string is attached. For a quark with infinite mass, r?, = 0; for 
finite mass = tq defined in equation (4.3). Therefore, we must distinguish 
the following two situations: 

► Infinite Quark Mass. Taking into account the definition (C.2) and the 
asymptotics (C.4) and (C.6), we impose: 

1 = v[/(0, u) = ^ hm A-2/3yl,^,„,,,, = (C.34) 



where in the last line we used the definition uog = uvg. 

Using the chain of equations (C.34), (C.33) and (C.16) we can finally fix, 

in the large u regime, the coefficient Ch as: 

\Ch\=(il{Wsr,fl^X%''uj. (C.35) 

From this one obtains the coefficient '^h defined in (4.9) using the relation 
^(x) = 7^-l/2^(x) and (C.8). This finally yields (4.39). 

► Finite Quark Mass. In an analogous way as for an infinitely massive 
quark, we use the form of the solution C.24 to write the normalization 
condition at the cutoff rq. However, the background scale factor h{x) 
cannot be approximated by the expression in (C.6) at a generic cutoff rg, 
but we rather need to keep and evaluate the full 7?."^/^(x) ~ b{x)/y/^ 
entering the expression for \&(x), without approximating it to \E'(x) ~ 
y/^xX{x)~'^^^. In fact, the necessary condition allowing to make use of 
this approximation is that A(x) ^ 1. The normalization condition then 
reads: 

1 = ^(rg, U) = TZ'^/^{Xq) [A^^lJsource{xQ) + A2^|J^ev{x q)] . (C.36) 

We now have to use the general form of the functions Fi and -F2, and keep 
both the source and vev solutions, since at rq the solution ip^ev in not 
negligible. Hence, using the relations (C.16) and (C.33) yielding 

A2 = iXti^Ai, 

the normalization condition determines Ai as 

Ai = -fUsXqTZ^^'^ {xq)[Fi{xQ) {cos{-fuJsx) + 7w^x sin(7a;^x)) 

+iXtl,^F2{xq) {sin{'-fUJsx) - '^Usxcosi'yujsx))] \ (C.37) 

Moreover the asymptotics of Fi and F2 must be generalized w.r.t. (C.25), 
since we need to evaluate them at xq , which is not necessarily close enough 
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to the boundary to allow us to use (C.25) — namely if A(xq) is not very 
small. In fact, one should keep the exact form of the solution at zero 
frequency (C.23), which reads 

C rl'r' 

i^hvi^) - C,V^n'/\x) + ^7^V2(a:) y ^ . (C.38) 

Comparing it to (C.24) one obtains the following behavior of function Fi 
and F2: 

(C.39) 

Consequently, one gets the value of the modulus square of the coefficient 



^8/3 -^2(^Q 



1 + {'JUsXqY 



+ ( ^ ^ ) (^HI^sXq) - 7W,XqCOs(7W,Xq))' 



(C.40) 



As we explained previously in this appendix, Fi and F2 are approxi- 
mated by (C.39) and (C.30), respectively for jUsX <^ 1 and x Xtp — 
a/2/7(X'. Therefore, substituting this asymptotics in \Ch\'^ and using equa- 
tion (4.17), we arrive at the results of section 4.4. 

Throughout the calculation we assumed that one can neglect the 0{X) terms 
compared to terms of order 0{1). This criterion is indeed satisfied in the numerical 
examples we study in this paper. In the cases where this is not satisfied, or one 
needs better accuracy in the WKB approximation, than one should work out the 
sub-leading corrections. 



D. Correlators in = 4 

In this Appendix we would like to collect and derive some results on the imaginary 
part of the retarded correlator for the Af = 4 theory. Some features of this quantity 
were discussed in [14] (see also [18]). More specifically, in [14] the symmetric corre- 
lator for J\f = 4 — related to Im Gr by equation (4.18) — is numerically computed 
and an analytic approximation is proposed. Here we show the numeric result for 
Im Gf{ and compare it to a linear plus cubic function. The advantage of considering 
Im Gr with respect to Gsym is the possibility of distinguishing the corrections to the 
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large and small frequency behavior associated to the coth factor in (4.18) from the 
wave function corrections, appearing in the coefficient in (4-17). In Ggym both 
kinds of corrections arise, while in Im Gr only the wave function contributes. 

Unlike the non conformal case, the = 4 correlators for the longitudinal and 
transverse modes only differ by a factor 7^, due to the fact that the wave functions 
satisfy the same equation for both kinds of modes. The extra 7^ comes from the 
factor in equation (4.1) for the longitudinal modes, that is constant in the 
conformal limit: Z — )■ I/7. Moreover, the fluctuation equation for both transverse 
and longitudinal modes only depends on the dimensionless variables x = r/r^ and 



UJ 



dx 



A 



x'^ ^ ^ , „ , (4a) 



-dj^{x,u}) 



~A2 



X2(l -X^) 

Following the same steps as in Section 4, we obtain: 

^R{i,Q) = Ciuj)il-xr'' + --- (D.2) 
2 



Im Gj, = -v/A^V7(7rT)=' cI;|C((I;)|2 (D.3) 

TT 

The WKB approximation implies that the imaginary part of the retarded cor- 
relator grows as for large frequencies, in the infinite mass case. More precisely, 
taking the conformal limit, b{r) £/r and i/ig ^^^=4, we obtain 

C{u) = 4u + 0{1) (D.4) 

and 

Im Gji = -f^Hm gI - — VAAr=4W^ for a; > -. (D.5) 

27r Ts 

This result^® can be obtained by applying the WKB method of the previous subsec- 
tion the the wave functions obeying equation (D.l). 

On the other hand, the result for small frequencies is well known [12, 14, 18] since 
it provides the diffusion constants^^. It is derived by analyzing the wave function 
solution to (D.l), in the regime where uvg = 4a) is small. The expressions for small 
frequencies read 

G{u) = l + 0{u;) (D.6) 



-"^^The fact that the large-w hmit is a cubic power-law can be expected from the corresponding zero- 
temperature result, in which one can analytically compute this quantity, and from the consideration 
that the large-frequency limit should be conformal. 

"'^^In reference [18] it is claimed that the imaginary part of the correlator is exactly linear in w, 
for all frequencies. This is not so, as explicitly shown here. 
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and 

Im Gji = j-Hm G{ ~ ^^Xj^=^T^uj, for w < -. (D.7) 

This result can also follow from taking the conformal limit of the diffusion constants 
derived in the next section. In fact, equations (5.7)-(5.8) are related to (D.7) by the 
formula (5.1). 

In figure 15 we show the numeric result for the TV = 4 correlator for the transverse 
modes compared to a linear plus cubic polynomial approximation, 

Im Gji ^ ciu + csw^ (D.8) 



with Ci and C3 given by the small and large frequency asymptotics of the correlator, 
as in equations (D.7) and (D.5), respectively: 




Figure 15: The picture on the left shows three curves (for \j\f=i = 5.5 and T = 250 MeV): 
i) the numeric result for Im G/j (red), ii) the linear plus cubic function cioj + c-iuj^ where we 
use (D.9) for ci and C3 (black dashed), iii) the difference of the linear plus cubic function 
with respect to the numeric result for Im G/j, multiplied by a factor of 20 (blue dotted). 
The right plot shows the relative difference between the linear plus cubic function and the 
numeric result for the correlator. 

The relative difference between the polynomial (D.8) and the numeric result 
for the imaginary part of the retarded correlator vanishes, as expected, both for 
small and for large frequencies. Nevertheless, the plots of figure 15 show that for 
1 < ojTs < 4 there is a sensible difference, of the order of 10-30%, between the two 
results^°. 

^°In Figure 15 we had to fix a value for T and \js/=i in order to plot Im Gr, but it is important to 
stress that the dependence on these quantities is simply given by the overall prefactor appearing in 
equation (D.3): the non-trivial part of the correlator depends on T only through the combination 
uj = Lu/Ts, and is independent of Aaa=4- 
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In the non-conformal cases the w — )■ oo hmit of the correlator is again controlled 
by the UV. Conformal invariance again implies an 0{u'^) behavior although it may 
be corrected by logarithms. 

It is also interesting to study the sub-leading corrections to IttiGr in the freqency. 
We first consider small u. The leading term is determined by the Kubo's formula to 
be linear. On the other hand in any P-invariant quantum field theory the imaginary 
part of a retarded correlator is guaranteed to be odd in u, see for example [57] 
for a recent discussion. Therefore we learn that the sub-leading correction at small 
frequencies is cubic. 

At high frequencies, the question is answered by extending the WKB analysis 
of Appendix C to sub-leading frequencies. First of all we note that the form of 
the Schrodinger potential near horizon as written in (C.9). Now, the WKB wave- 
functions are obtained from this by making use of (C.12) and (C.13). The piece that 
is relevant for the current discussion is the l/^/p part that in front of the cosine 
and the sine terms that clearly yields a 0{uj~'^ correction to the wave-functions in 
the large u limit. Matching the WKB solution to the near-horizon solution as in 
Appendix C, yields C2 = iCi and one clearly obtains, 

Ch(xCi{l + 0{u~^)). (D.IO) 

as a correction to equation (C.16). 

On the other hand, the Schrodinger equation near the boundary, in the TV = 4 
theory, can be written as a function of xUg- Therefore the wave-function near the 
boundary is a function of xUg- This means that, upon demanding unit normalization 
of the source term — the coefficient Ai in Appendix C — one always has x l/o; ~ 1 
with no correction perturbative in Therefore matching the wave- functions near 
the boundary and the WKB region yields Ci oc oc u. Substituting in (D.IO) we 
learn that ~ + 0{u~^). Finally, using (4.17) that determines the frequency 
dependence as ImGn oc a;|C;ip, we find that the subleading correction is linear in u, 
in the high u limit. 
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